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INTRODUCTION 

The promise of filamentary composite materials, whose 
development may be considered as entering its second genera- 
tion, continues to generate intense interest. Such interest 
is well founded, since it is based on the possibility of 
using brittle materials with high modulus, high strength, 
but low density in composites which fail in a non-catastroph- 
ic manner. Such fiber reinforced composite materials offer 
substantially improved performance and potentially lower 
costs for aerospace hardv/are. 

Much progress has been achieved since the initial de- 
velopments in the mid 1960 's. Only a very few applications 
to primary aircraft structure have been made, however, and 
those are in a material-substitution mode and - with the 
exception of experiments only just now being conducted on 
large passenger airplanes - on military aircraft. 

To fulfill the promise of composite materials more com- 
pletely, requires a strong technology base. NASA and AFOSR 
have realized that to fully exploit composites in sophisti- 
cated aerospace structures the technology base must be im- 
proved. This, in turn, calls for expanding fundamental 
knowledge and the means by which it can be successfully ap- 
plied in design and manufacture. Not the least of this 
effort is to learn how to design structures specifically to 
capitalize on the unique properties of composite materials. 
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It also calls for expanding the body of engineers and scien- 
tists competent in these areas. As part of their approach 
to accomplishing this, NASA and AFOSR have funded the current 
composites program at Rensselaer. The purpose of the RPI 
composites program is to develop advanced technology in the 
areas of physical properties, structural concepts and analy- 
sis, manufacturing, reliability and life prediction. Con- 
commitant goals are to educate engineers to design and use 
composite materials as normal or conventional materials. A 
multifaceted program has been instututed to achieve these 
objectives. The major elements of the program are: 

1. CAP COMP (Composite Aircraft P^rogram Comp onent) . 
CAPCOMP is primarily a graduate level project being conducted 
in parallel with a composite structures program sponsored by 
NASA and performed by a private, aerospace manufacturing con- 
tractor. The first component redesign is being done in con- 
junction with the Boeing Commercial Airplane Company. The 
main spar/rib region on the Boeing 727 elevator, near its 
actuator attachment point, was selected, with Boeing's advice 
and the concurrence of NASA/AFOSR, for study in CAPCOMP. The 
magnitude of the project - studying, designing, fabricating 
and testing of a relatively small but highly stressed region 
on the elevator - is both consistent with Rensselaer's capa- 
bilities and a significant challenge. The selection of a 
portion of a full-scale flight hardware structure assures 
relevance to this project's direction. 
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Visits to Boeing were conducted by Professor Hoff and 
several of his students, and the first serious design work 
began on two alternative designs. Each was pursued to the 
point of preliminary analysis and testing. One of these was 
selected for more detailed analysis, redesign, complete fab- 
rication and testing. Completion of fabrication facilities, 
the first complete part and test fixture development are re- 
ported in Part I. 

2. CAPGLIDE (Composite Aircraft Program Glider ) . This 
undergraduate demonstration project has as its objectives 
the design, fabrication and testing of a foot-launched, 
ultralight glider using composite structures. A flight ve- 
hicle was selected to maximize student interest and to pro- 
vide the students with a broad -based engineering experience. 
For those students continuing with graduate work at PJ>I , 
CAPGLIDE is intended to provide natural progression to CAP- 
COMP. The progress on the CAPGLIDE project to date has been 
very good. Seven professors and approximately 30 students 
were actively engaged in the project during this period; 
that is, the Spring '80 semester and summer period. High 
point of the project to date was initial flight testing of 
the glider, dubbed "RP-1", and first public flight demon- 
stration on September 16, 1980. A description of the status 
of the CAPGLIDE project at the end of the current reporting 
period is given in Part II. 

3. COMP AD ( Compu ter Aided Design) . A major thrust of 
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the composites program is to develop effective and efficient 
tools for the analysis and design of composite structures. 
Rensselaer and NASA Langley have jointly implemented the use 
of the SPAR code on minicomputers, and the work at Rensselaer 
has made "virtual memory" available to those using SPAR. 
Special RPI software has been developed to make existing 
programs more useful, particularly as regards graphics. At- 
tention for the current period has focused on preprocessor 
and postprocessor developments; details are reported in Part 
III. 

4. Composites Fabrication and Test Facility . Struc- 
tural design engineers, educated only by course work and 
design projects limited to paper, often fail to sense or 
appreciate problems involved in fabrication. The actual 
fabrication and testing of composite structural components 
provides this training and the final validation for the de- 
signs in our CAP projects. RPI's Composites Fabrication and 
Test Facility is located in the laboratory and high bay areas 
of the Jonsson Engineering Center. Equipment is available 
for compression molding parts as large as 19" x 19" and vac- 
uum bagging parts of much larger size. Panels, approximately 
4' X 20', have been made by vacuum bagging. A pressure ves- 
sel for long narrow parts, such as spars, has been designed 
and built. All of these techniques and facilities have been 
used in the CAP GLIDE part of the project, described in Part 
II. A second vessel capable of fabricating shorter but wider 
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parts at higher pressures was designed and completed in the 
current period. More information on this development is con- 
tained in Part I, CAPCOMP. 

5. INSURE (Innovative and Sup porting Research) . The 

criteria for selecting research projects to be conducted 
under this program include the following: (a) they must 

anticipate critical problem areas which may occur in the CAP 
or NASA/AFOSR programs or (b) solutions to existing problems 
are not yet satisfactorily in hand. During the reporting 
period six such projects were conducted as part of the 
program. Results from the ongoing projects are reported in 
Part IV. 

6. Curriculum Revisions . The goal of educating engi- 
neers to think of composites as normal or conventional mate- 
rials has required changes in curriculum. Since the initi- 
ation of this program, almost all Rensselaer engineers take 
introductory courses which incorporate the concepts of ani- 
sotropy and composite materials. In addition, six special- 
ized courses in composites have been offered during the past 
three years to develop those special skills required of 
students involved in the composites program. A new course 
was introduced in the Fall of '78 semester on composite de- 
sign and analysis using programmable hand calculators, a 
central mini and full frame computers. A new graduate level 
advanced topics course with the title "Advanced Finite Ele- 
ments" was offered for the first time in September 1979. 
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The additions of the SAP and SPAR computer codes and 
the growing availability of interactive computer graphics 
under our COMPAD program element have reached the point 
where our engineering students are using these facilities as 
everyday working tools for design, analysis and visualization 
purposes. We have thus achieved one of the principal goals 
of the curriculum development activities. 

Two continuing education special courses, outgrowths of 
the composite materials and structures program, were pre- 
sented during the month of July. They were "Composite Ma- 
terials & Structures" and "The Finite Element Approach in 
Structural Analysis". Each course lasted a week, and they 
were run serially so that a participant might come for both 
on the same trip. The level of the material was planned to 
be suitable for graduate engineers and useful to managers of 
engineering structures activities who are involved in tech- 
nical work but who may not have taken courses for several 
years. Because of the v/ide variety of special courses 
available throughout the United States dealing with this 
subject matter, only rather unique aspects justify additional 
offerings in these areas. Our programs were planned to be 
unusual in the respect that "hands-on" experiences were in- 
herent and required for completion of each course. Use of 
hand-held programmable calculators, application of optimiza- 
tion programs using computer graphics, lay-up and cure of a 
simple part in graphite epoxy and testing tensile, shear and 
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compression coupons to failure were all part of the first 
course. As part of the second course, workshop sessions 
were held in KPI’s Computer Graphics Center. In these ses- 
sions students defined, "constructed" and analyzed their own 
finite element models on refresh graphics terminals. 

Faculty-staff and student attendees for the courses are 
as follows; 


Composite Materials & Structures 

Program Director: 

R. G. Loewy 

Institute Professor 

Rensselaer Polytechnic Institute 


Faculty: 

R. Judd Dief endorf 
Professor, Materials Eng. 
Rensselaer Polytechnic Inst. 

H. Gunter Helwig 
Asst. Research Professor 
Aeronautical Eng., Mechanical 
Eng. and Mechanics 
Rensselaer Polytechnic Inst. 

B. Walter Rosen 
President 

Materials Sciences Corp. 
Students: 

Cecil G. Armstrong 
Research Engineer 
Rensselaer Polytechnic Inst. 

David J. Brangwin 
Senior Engineer 
Harris Corporation 

Richard H. Cole 
Mechanical Engineer 
Benet Laboratory 
Watervliet Arsenal 


Stephen VI. Tsai 
Chief, Mechanics and Sur- 
face Interactions Branch 
Air Force Materials Lab. 

Wr ight-Patterson AFB 

Dick J. Wilkins 

Sr. Engineering Specialist 

General Dynamics 


Michael J. Dennis 
Specialist II, Structures 
Technology 
Boeing Vertol 

John H. Kelbel 
Mechanical Engineer 
Nurad Inc . 

William L. Poesch 
Member of Technical Staff 
MRJ, Inc. 
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J. B. Senior 

Project Stress Engineer 

Rolls-Royce Ltd. 

Derby, UK 


Alan Servoz 
Engineer 
Aerospatiale 
Toulouse, France 


The Finite Element Approach in Structural Analysis 


Program Director: 

Mark S . Shephard 
Assistant Professor 
Civil Engineering 
Rensselaer Polytechnic Inst. 

Faculty: 

Richard H. Gallagher 

Dean, College of Engineering 

University of Arizona 

John A. Swanson 
President 

Swanson Analysis Systems, Inc. 

Students : 

Richard Adkins 
Design Engineer 
Ashland Petroleum Co . 

Frank P. Hill 
Manager, Analytical Eng. 
Combustion Engineering Inc. 

John P. Housprup 
Instructional Consultant 
Combustion Engineering Inc . 

Robert Jankowski 
Sr. Structural Engineer 
Brumman Aerospace 

Edward Kirchner 
Research Engineer 
Farr el Corporation/Div. USM 

William Lewis 
Mechanical Engineer 
Mechanical Technology Inc . 

Sudhir C. Mehrotra 
Research Scientist 
Vigyan Research Associates, 

Inc , 


Gary S. Miller, Sr. 

Staff Engineer 

General Physics Corporation 

Joyceann Montini 
Associate Engineer 
Boeing Vertol 

Charles A. Pagel 

Mechanical Engineer 

Naval Weapons Support Center 

Alan Servoz 
Engineer 
Aerospatiale 
Toulouse, France 

John St imp son 
Mechanical Engineer 
Raytheon Company 

John Talbourdet 
Structural Engineer 
Honeywell Electro Optics 
Systems Center 

Roger Weseloh 

Sr. Design Engineer 

Westinghouse 
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Reaction to these offerings on the part of the attendees 
was highly favorable, and we plan to give the courses again 
in 1981. 

7 . Technical Interchange . 

a) Mr. Hiromitsu Otsuka, who joined the faculty at RPI as a 
Visiting Research Associate in Mechanical Engineering, Aero- 
nautical Engineering and Mechanics, returned to Japan in July 
1980. Mr. Otsuka is a member of the staff of the Industrial 
Products Research Institute of the Ministry of International 
Trade and Industry, where he was engaged in developing a 
light weight wheelchair of fiber reinforced plastic materials. 
At RPI he conducted acoustic emission studies of composites 
under load, together with RPI's Professor H. Scarton, during 
this period. The results of this work are reported in Part 
IV, INSURE. 

b) Technical Meetings: Technical meetings, on- and off- 

campus, provide important opportunities for interchange of 
technical information. Because of the large number of com- 
posites meetings, a central catalog with all upcoming meet- 
ings is being maintained and distributed periodically. In 
this way we help assure that a Rensselaer staff member will 
participate in important meetings. The calendar for this 
reporting period is shown in Table 1. Meetings attended by 
RPI composites program faculty/staff during the reporting 
period are shown in Table 2. Some meetings particularly 
relevant to composites, held on-campus with off-campus 
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TABLE 1 

CALENDAR OF COMPOSITES -RELATED MEETINGS 
for the period May '80 to September '80 


5/6-8 

5/12-14 

5/13-15 

5/14-15 

5/15-16 

5/19-23 

5/27-30 

5/28-29 

6/2-5 

6/4-8 

6/30-7/4 

8/11-13 

8/17-23 

8/25-30 

9/1-5 

9/2-3 

9/3-5 

9/10-12 


25th National SAMPE Symposium and Exhibition, San 
Diego, CA. "Sponsored by SAMPE." 

21st Structures, Structural Dynamics and Materials 
Conference, Seattle, WA. "Sponsored by AIAA/ASME." 

Annual Forum, Washington, D. C. "Sponsored by AHS." 

Impact for the '80's (Materials/Structures Panel), 
Cleveland, OH. 

Analysis of USAF Aircraft Engine Structural Dura- 
bility and Damage Tolerance, Seattle, WA. "Sponsor- 
ed by AIAA." 

4th International Conference on Pressure Vessel 
Technology, London, England. "Sponsored by ASME." 

lUTAM Symposium on Physical Non-Linearities in Struc- 
tural Analysis, Senlis, France. 

M.A.B. Meeting: Committee on High Temperature Metal 

and Ceramic Matrix Composites, Washington, D. C. 

lUTAM Symposium on Three Dimensional Constitutive 
Relationship and Ductile Fracture, Dourdan, France. 

198 0 SIALl National Meeting, Alexandria, VA. 

Carbon '80, Baden-Baden, Germany. 

7th Atmospheric Flight Mechanics Conference (includes 
Aeroelastic Tailoring of Flexible Aircraft), Danbers, 
MA. "Sponsored by Battelle, Columbus Labs." 

15th International Congress of Theoretical and Ap- 
plied Mechanics, Toronto, Canada. "Sponsored by AIAA." 

3rd International Conference on Composite Materials, 
Paris, France. "Sponsored by AIAA." 

8th International Conference on Rheology, Naples, 

Italy. 

International Conference on Vibrations in Rotating 
Machinery, Cambridge, England. "Sponsored by ASME." 

USA-USSR Symposium on Fracture of Composite Materials 
Bethlehem, PA. 

Mechanics of Non-Destructuve Testing, Blacksburg, VA. 
"Sponsored by ARO-D; NSF and V.P.I." 
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TABLE 1 (continued ) 


9/15-17 

9/15-19 

9/16-17 

9/16-17 

9/17-19 

9/24-27 

9/30 


1st Chautauqua on Finite Element Modeling, Harwich- 
port, MA. 

Composite Materials Processing and Quality Assur- 
ance, Los Angeles, CA. 

4th Meeting, Commission on Acoustic Emission from 
Reinforced Plastics, Montery, CA. "Sponsored by 
Lawrence Livermore National Lab. and Acoustic Emis- 
sion Tech. Corp." 

6th European Rotorcraft and Powered Lift Aircraft 
For\im, Bristol, England. 

IPAD National Symposium, Denver, CO. "Sponsored by 
NASA." 

4th International Symposium on Composite Materials, 
Rydzyna, Poland. 

Carbon/Carbon : Composite Technology, Dallas, TX. 


9/30-10/2 Army Symposium on Solid Mechanics, Yarmouth, MA. 

"Sponsored by Army Materials and Mechanics Research 
Center . " 
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TABLE 2 

COMPOSITES -RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
for the period May '80 to September '80 


5/11-15 

5/27-30 

5/28-29 

6/2-5 


6/4-8 

6/30-7/4 

8/25-30 

9/1-5 

9/15-17 

9/16-17 

9/24-27 

9/30 


21st Structures, Structural Dynamics and Materials 
Conference (Prof. Helwig) , Seattle, WA. 

lUTAM Symposium on Physical Non-Linearities in 
Structural Analysis (Prof. Diefendorf ) , Senlis, 
France. 

M.A.B. Meeting: Committee on High Temperature Metal 

and Ceramic Matrix Composites, (Prof. Diefendorf) , 
Washington, D.C. 

lUTAM Symposium on Three Dimensional Constitutive 
Relationship and Ductile Fracture (Prof. Krempl) , 
Dourdan, France. 

* 

1980 SIAM National Meeting, (Prof. Brunelle ), Alex- 
andria, VA. 

Carbon '80 (Prof. Diefendorf), Baden-Baden, Germany. 

3rd International Conference on Composite Materials 
(Prof. Diefendorf) , Paris, France. 

8th International Conference on Rheology (Prof. 
Diefendorf) , Naples, Italy. 

1st Chautauqua on Finite Element Modeling (Profs. 
Feeser and Shephard**) , Harwichport, MA. 

4th Meeting, Commission on Acoustic Emission from 
Reinforced Plastics (Prof. Scarton) , Montery, CA. 

4th International Symposium on Composite Materials 
(Prof. Diefendorf), Rydzyna, Poland. 

Carbon/Carbon: Composite Technology (Prof. Diefen- 

dorf) , Dallas, TX. 


Presented a paper titled, "The Use of a Double Affine Trans- 
formation to Obtain a Class of Generalized Elastodynamic So- 
lutions for Orthotropic Plates" . 

Presented a paper titled, "RPIFEP - An Interactive Finite 
Element Package" . 
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speakers, are listed in Table 3. A list of composite-related 
visits to relevant organizations by RPI faculty/staff/stu- 
dents, with the purpose of each visit outlined, is presented 
in Table 4 . 

In summary, the NASA/AFOSR Composites Aircraft Program 
is a multi-faceted program whereby aeronautical, mechanical 
and materials engineers must interact to achieve its goals. 
"Hard-nosed" engineering of composite aircraft structures is 
balanced against research aimed at solving present and future 
problems. In the following sections, detailed descriptions 
of the CAPCOMP, CAPGLIDE, COMPAD and INSURE programs are pre- 


sented . 
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TABLE 3 

COMPOSITES -RELATED MEETINGS /TALKS HELD AT RPI 
(May '80 - September '80) 


l2Ei£ Date 

F-16 and Other Developments 7/17/80 
at General Dynamics 


Bending-Torsion Coupling 7/29/80 

Effects in Cascade Flutter 


A Hierarchy of Conforming 9/22/80 

Timoshenko Beam Elements 

Oscillating Cascade Fluid 9/30/80 

Energy Extractor 


Speaker (s) 

Dick J. Wilkins 
Sr. Engin. Specialist 
General Dynamics 
Fort Worth, TX 

Oddvar O. Bendiksen 
Dir. of Project Engin. 
Pacific Airmotive Corp. 
Burbank, CA 

Alexander Tessler 
Northrop Aircraft 
Los Angeles, CA 

Lee Arnold 

Emeritus Professor of 
Applied Science 
New York University 
New York, NY 


TABLE 4 

COMPOS ITES-P^ELATED VISITS TO RELEVANT ORGANIZATIONS 
by RPI Faculty/Staf f/Students 


Visited 


Date by Prof (s) 


Purpose 


G. E. Space Di- 9/25/80 

vision 

Dr. Carl Zweben 


R. G. Loewy To review composite 
structures activ- 
ities 


Drexel University 9/25/80 R. G. Loewy 

Professors: 

P . C . Chou 
D. Mortimer 
J. Awerback 
A. Wang 
J. Rose 
G. Harlow 


To discuss programs 
in composites re- 
search and visit 
laboratories 
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CAPCOMP (Gomposite Aircraft Program Com.ponent) 

(N. Hoff, R. Loewy) 

CAPCO^^P is a program to design flight critical struc- 
tures to take the maximum advantage of composite materials. 
By combining the efforts of experienced faculty with bright 
and well trained but inexperienced graduate students in an 
environment relatively free of traditional design and manu- 
facturing processes, we intend to devise new and hopefully 
useful design concepts. 

There is sufficient information available to prove that 
many structural elements can be made lighter by using ad- 
vanced composites today than by using metals. But if such 
elements have to be joined by methods other than adhesive 
bonding, difficulties and uncertainties arise which can be 
eliminated only through conservative designs with their at- 
tendent penalties in weight or by extensive and expensive 
programs of "cut and try" . This stands as one important im- 
pediment to full adoption of composites by the aerospace in- 
dustry. 

On the basis of these considerations Rensselaer Poly- 
technic Institute began, as the first task aimed at new 
structural concepts, the design using composites of a joint 
used in an airplane elevator. To make the design realistic, 
an existing metal airframe component was chosen for redesign 
in composites. The existing design chosen was that of the 
Boeing 727 elevator actuator attachment. We conceive of 
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this work as carrying forward a Structures Demonstration 
Program using the actuator attachment joint of the 727 ele- 
vator, paralleling that of NASA and its aerospace contractor, 
the Boeing Commercial Airplane Company. Our design, fabri- 
cation and test effort emphasizes designs using advanced 
composite construction for the purpose of minimizing the 
weight of the structure, but on a scale consistent with the 
university context and funding level. The staff of RPI is 
very grateful to the Boeing Company and its engineers for 
their wholehearted support of this work at RPI. 

During the last reporting period, one of two different 
designs reported previously (July and December 1979) as 
suitable for replacing the largely metal attachment produced 
by Boeing was selected for further analysis, redesign, fabri- 
cation and testing. We call it the "Berg Design", and it 
makes use of quasi-isotropic , graphite-epoxy laminates. The 
second design, called the "Muser Design", which made a de- 
liberate attempt to use uniaxial graphite-epoxy tape to as 
great an extent as possible, was suspended. It has led, how- 
ever to more generic research efforts to maximize the load 
carrying ability of pin-loaded holes in composite membranes 
and plates. These efforts are directed at improving the de- 
sign of the attachment of the rib flanges to the elevator 
skin through mechanical (i.e., pin-type) fasteners and at 
maximizing the efficiency of the movable, heavily loaded, 
actuator attachment points in CAPCOMP. In addition, however. 
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and increasingly, these generic effects are motivated by the 
general and widespread utility that will result from improve- 
ments in pin-loaded hole designs for composites. 

These efforts have been jointly directed by Dr. Nicholas 
J. Hoff, in part time arrangement and by Dr. Robert G. Loewy. 
They are further described in the following paragraphs. 

1. The Elevator and Its Attachment 

To place the progress reported here in proper context, 
some descriptions used in earlier reports are repeated. The 
conventional aluminiim alloy elevator of the Boeing 727 is 
shown in the upper half of Figure I-l. The lower half of 
the figure is the new version of the elevator redesigned by 
Boeing in graphite epoxy; it is evident from the pictures 
that the latter is composed of fewer parts than the former. 
However, the actuator fitting of the new design is still 
manufactured of aluminum alloy. This fitting is shown in 
Figure 1-2. The fitting is attached to outboard and inboard 
portions of a new graphite-epoxy spar and to a graphite- 
epoxy nomex-honeycomb rib as indicated in Figure 1-3. 

The attachment was designed by Boeing to carry loads up 
to 19,000 pounds. The direction of the load varies as the 
elevator rotates over an angle of 28 degrees from the full- 
down to the full-up position. 

2 . Berg Vs. Design 
(R. Loewy) 

Berg's design for the elevator rib is shov/n in Figure 1-4 





BOEING ELEVATOR ASSEMBLY 


Figure T-1 









COMPOSITE ELEVATOR ACTUATOR ASSEMBLY (Boeing Design) 

Figure 1-3 




ACTUATOR LUG TO RIB TRANSITION - TOP VIEW 


Figure i-4 
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and 1-5. The first of these figures compares the Boeing com- 
posite design with Berg's design; both are shown. Figure 1-5 
shows both how the spar and new actuator rib would be assem- 
bled and also the build-up of thickness in the lug areas an- 
ticipated as needed to provide the required bearing strength. 
Note that the edges of some of the layers are bent 90 degrees 
to form flanges to which the upper and lower skins of the 
elevator can be attached. Attachment of flanges to elevator 
skins will be by means of titanium Hi-Loc fasteners. The 
right-hand and left-hand graphite-epoxy webs are stabilized 
by a layer of nomex-honeycomb between them. 

The design analyses described in the previous reports 
led to an elevator attachment rib with a variable thickness 
web. The distribution of thickness was chosen using succes- 
sive finite element calculations to produce a relatively 
uniform stress level. It is shown in Figure 1-6. 

During the current period a new pressure vessel, design- 
ed and fabricated especially for fabricating carbon-epoxy 
components in the CAPCOMP size category, was completed. Its 
useful dimensions are 7 '8" x 18" diameter, gage pressures of 
90 psi are available from the compressed air supplies in our 
Composites Fabrication Laboratory, and internal temperatures 
of up to 140°F can be generated by electrical conduction 
heating elements powered by a Thermae cure-cycle controller. 
Figure 1-7 is a photograph of the completed pressure vessel. 

The Berg Design, composite 727 aircraft elevator actuator 



PROPOSED SPAR AND RIB ASSEMBLY (Berg Design) 


Actuator Rib 


Elevator Spar' 


Elevator Hinge 







NOTE; Two luge on actuator rib 

are inserted through slots 
in web of elevator spar 
for assembly. 


Figure 1-5 
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Figure 1-6 









gglii 


1 







30 


attachment rib is shown in various stages of assembly tsee 
Figures I-8a, -8b and -8c) . 

As noted in the last report, the Boeing Commercial Air- 
plane Company provided a full description of the testing 
procedures they used to qualify their 727 elevator actuator 
attachment structure. Figure 1-9 shows the Boeing test set- 
up; note that two hydraulic cylinders (mounted horizontally) 
provide the actuator load and reaction through a vertical 
bar which is bolted to the two elevator actuator attachment 
points. In the absence of such mobile hydraulic load gener- 
ators at RPI, we plan to generate the equivalent loadings in 
an Instron testing machine through linkages as shown in Fig- 
ure I-IO. In the campus test set-up, loads will be applied 
vertically to a bar oriented horizontally which transmits 
the loads through two bolts to the elevator actuator attach- 
ment lugs. This linkage arrangement, which provides the 
equivalent of the Boeing load test, has been stress-checked 
and fabrication has been completed during this reporting 
period. 

Also shown in Figure I-IO is the "strong-back" structure 
(items 9, 10, and 11) and the reinforcements of the elevator 
rib (items 12 and 13) to ensure reaction-load transfer to 
the "strong-back" structure. These test fixtures (i.e., 
items 9, 12 and 13) and the aft elevator skin reinforcements 
(item 14) are essentially those used in the Boeing test. 

Their fabrication also was completed in this current report- 
ing period. 










Figure I-8a 


Left-Hand Member Elevator Actuator Attachment Rib 
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Elevator Attachment Rib Assembly Counted and 
Instrumented for Testing 
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From equilibrium of moments in item 5 in Figure I-IO 

(see free-body diagrams in Figure I-ll) , it is clear that 

P’Jl = [P/2]3£. Thus P = 2/3P*. To achieve full design load 

# 

on the elevator actuator rib attachment point (19,000 ), 

# 

therefore, the RPI Instron (maximum capacity 28,000 ) need 
develope only P = 2/3(19,000) =12,666^. Further, as de- 
signed, the moment equilibrium on the "strong-back structure" 
(items 8, 9, 10, 12 and 13) will be achieved, since this re- 
quires that M = P'fi, = 3/2P5, = PA, and, in fact, 3/2£ = 6.638" 
and A = 6.638". 

We expect the first proof load tests to failure of the 
Berg design actuator attachment rib to be conducted early in 
the next reporting period. Both strain gage and acoustic 
emission data will be taken during these tests. 

3. Optimizing Fiber Orientations in the Vicinity 
of Heavily Loaded Joints 

(N. Hoff) 

A. Introduction 

The purpose of this work is to find ways of arranging 
fibers in a composite structure that lead to high load- 
carrying capacity-to-weight ratios in the neighborhood of dis 
continuities such as mechanical joints and openings. The 
paucity of theoretical information of this kind and the scar- 
city of long-time practical experience with composite joints 
have compelled designers to use metals for fittings rather 
than composites in most highly loaded cases. It is believed 
that composites offer the possibility of saving weight in 
such joints and thus merit a detailed investigation. 
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B. Work Completed before May 1 ' 

As described in the previous Semi-Annual Progress Re- 
port, a manufacturing technique has been developed for tensile 
specimens incorporating a hole, with fiber arrangements vary- 
ing with distance from the hole in a manner intended to re- 
duce stress concentration factors. A few specimens were built 
and tested to failure. 

A start was made toward the establishment of a theory of 
the stress distribution around empty holes when uniform uni- 
axial tension is applied to the specimen far away from the 
hole. Such a theory has existed for many years in the case 
of orthotropy in a Cartesian system of coordinates, but no 
solutions have been known in the case when the plate is ortho- 
tropic in a polar coordinate system (cylindrical orthotropy) . 
Yet this kind of orthotropy promises to reduce the stress con- 
centrations. 

Our earlier investigations discovered the fact that a 
uniform state of tension cannot exist, even far av/ay from 
the hole, in a cylindrically orthotropic plate under edge 
loading when properties are independent of distance from the 
hole, except in the special case when the compliances are 
connected by the equation 


- “re ' ^66 ‘ ° 

where =r6 = ®re/®9e’ ®66 = ®66/®O0- 


(In the middle of pg. 54 the term Sgg was written as S 


00 


* 

Pg. 54 of 38th Semi-Annual Progress Report on Composite 
Structural Materials. 
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through a typing error.) In this case it was found possible 
to derive a closed-form solution for the circumferential 
stress at the edge of the hole. This solution is the next 
equation on pg. 54 {see footnote on previous page). No solu- 
tion of the problem was given in the general case. 

C. Progress Made during Reporting Period 

During the present reporting period our effort was con- 
centrated on obtaining analytical solutions for the case 
where properties are independent of distance from the hole, 
and we have made significant progress in this direction. 
Complete closed-form solutions have been found for the 
stresses around the hole for any arbitrary combination of the 
values of the compliances S , S and S,, as long as 

these values are constant. 

One interesting result of the analysis is that the four 
compliances can be combined into two parameters in such a 
manner that the stress concentration factor becomes a func- 
tion only of these parameters and not of the four individual 
compliances. One of the parameters is the ratio of the 

compliances for extension in the radial and circumferential 
directions: 

S = S /S.-. 
rr rr' 0 0 

The other parameter S is defined as 

+ =66 = 2(s^e/=ee> + 

The earlier conclusion that a uniform state of stress 
cannot exist in the cylindrically orthotropic plate, except 
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when the values of the compliances are related according to 
a particular condition, has been confirmed. 

The analytical work shows that the solution depends 
strongly on the values of the two parameters defined earlier. 
In the first range of the parameters the stress function and 
the expressions for stress and strain are exponential func- 
tions of the radial distance r with real exponents. Natur- 
ally these can be combined into hyperbolic functions of real 
argument. 

In another range of the parameters the exponents are 
complex and the solutions are products of hyperbolic func- 
tions with trigonometric functions of real argument. Final- 
ly, in a third range the solutions are products of trigono- 
metric functions of real argument . 

Even within these ranges the stress-concentration factors 
show different patterns of behavior. In one subrange the 
factors remain finite when the outer radius R of the circular 
plate analysed increases beyond all bounds. In a second 
subrange the factors approach zero signifying that the trac- 
tions applied to the outer edge of the plate are in equilib- 
rium in a boundary layer parallel to the loading and excluding 
the hole and cause no finite stresses in the vicinity of the 
central circular hole. 

In a third subrange the stress-concentration factor 
approaches infinity if the outer radius R of the plate is 
allowed to increase beyond all bounds. Although strange at 
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fx 3 Tst iTGading', tiixs bsiisvxoir xs not contsTSiry to any of tiiG 
laws of mechanics. After all, infinite plates are a mathe- 
matical fiction. In the physical world no one has yet built 
an infinite plate, nor has he loaded it at infinity. 

The manuscript of a full report on this phase of our 
investigations has been completed and is now ready for typing 
in its final form. It is entitled, "Stress Concentrations in 
Cylindrically Orthotropic Composite Plates with a Circular 
Hole". In addition to the mathematical derivations, the 
report also contains physical explanations of the various be- 
havior patterns observed. 

A good start has also been made with an investigation 
of the stress-concentration factors in cylindrically ortho- 
tropic plates when the compliances are functions of the ra- 
dial distance r. In this case the compliances have the form 

* * 2h 

S. . = S . . + S . . e.r 

xj X3 X3 X 


where h is an integer and and e^ are determined by the 

prescribed properties at the hole (r =1) and at the outer 
edge (r = R) . The solution for the stresses in such a plate 
is based on the procedure described earlier, where the com- 
pliances are not a function of the radius r, and it involves 
power series. To test the calculation method, the distribu- 
tion of circumferential stress along the radius was determined 
using a pocket calculator. Compliances were chosen as fol- 
lows: at the hole (r = 1) values are set for composites with 

±45° fibers only; at the outer boundary (r = R = 10) the 
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properties are made those of a nearly quasi-isotropic material 
v;ith a fiber arrangement (02^/45®^/-45°g/902g) • Since a state 
of plane stress is assumed the stacking sequence is not im- 
portant for the calculations. Two composite materials were 
chosen for the examples. One case is a typical glass epoxy, 
where the compliances do not change greatly from r = 1 to 
r = 10, and the other is a typical graphite epoxy where the 
compliances exhibit a large variation in the end points, that 
is the values at r = 1 as compared to r = 10. Values of com- 
pliance between the end points were not calculated based on 
variations of fiber lay-ups as might be expected in practice. 
Rather, simple power series were picked which could be ex- 
pected both to be tractable and, at the same time, provide 
for realistic values. Figure 1-12 shows for the property 
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this behavior of the compliances. Two different varia- 


tions of are shown for glass epoxy and one for graphite 

epoxy. All values are nondimensionalized . The overall ge- 
ometry of the loaded plate is shown in Figure 1-13. 

The distribution of circumferential stress along the 
radius r at 0 = 'it/ 2 for glass epoxy is plotted in Figure 1-14. 
The two curves for a different variation of S.. with the ra- 

±3 

dius r, may be compared with the one where the plate consists 
of ±45° fibers only. The difference in the stress distribu- 
tion occurs mainly in the vicinity of the hole and somewhat 
at the outer boundary. The stress concentration factors 
differ from the value of 3 encountered with isotropic mate- 
rials, as expected. 
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The curves for circumferential stress in the case of 
graphite epoxy are shown in Figure 1-15. Note that when the 
compliances are a function of the radius r , the stresses do 
not differ greatly from a uniform average stress distribu- 
tion. The two stress concentrations, here, are smaller than 
unity, and the highest stress occurs not at the hole but at 
a point one or two radii away from the hole. These maximum 
stress concentration values are about 1.5. It is noted that, 
in the case of ±45° fibers throughout, the stresses near the 
hole have a negative sign, indicating compression. 

A more thorough treatment of this problem v/here the com- 
pliances are a function of the radius will be the subject of 
follow-up research. We consider this study of great impor- 
tance because, as noted earlier, in the optimal plate most 
likely the fibers near the hole will be arranged differently 
from those far away from the hole. 

D. Plans for Upcoming Period 

1. The manufacturing effort and the testing of 
specimens described in the previous report 
will be continued. 

2. Graphs will be prepared to show the stress- 

concentration factors as a function of the 
two parameters and S. Values will be 

computed from the expressions already de- 
rived; they should prove useful in practical 
design. 

3. The investigation of the stress-concentra- 
tion factors will be continued for cases 
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when the compliances (and thus the para- 
meters and S) are functions of the ra- 

dial distance r. It is expected that a 
technical report covering this Work will 
become available by the end of the up- 
coming period. 

4 . Supporting Research on Lightly Loaded Mechanical Joints 

(R. Loev/y) 

A. Introduction 

The work of graduate student Wonsub Kim continued to 
concentrate on stress distributions around pin-loaded holes 
in composite membranes and plates with no special fiber ar- 
rangements around the hole and the associated bearing strength 
and failure mechanisms. 

B. Status 

Three different approaches outlined in the previous re- 
port are being pursued in parallel, as a way to increase 
understanding of the behavior of plates subjected to bearing 
loads. Of these, emphasis during the past months was directed 
to finite element numerical computations and correlation with 
photoelastic analysis. 

It will be recalled that the finite element model shown 
in Figure 1-16 has been postulated. A constant strain, tri- 
angular shaped element is used and a 345 element model with 
212 nodes was generated by the program. The plan we formu- 
lated was to use several constraint equations to simulate the 
loaded edge condition. If we can match the strain distribu- 


tion from one of the finite element solutions with the 




Figure 1-16 Finite' Element Idealization and Boundary Conditions 
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fringes from the photoelastic experiment, we will then Icnow 
the contact zone assxamed between the loading pin and the 
edge of the hole in the corresponding computation was cor- 
rect. If successful, this will be a linear solution to an 
essentially non-linear problem. 


C. Progress During Reporting Period 

To put our model for the finite element analysis of a 
pin-loaded composite plate into the main analysis program, 
considerable time has been spent in data preparation. The 
output and input of the main program are in terms of element 
coordinates. It was, therefore, necessary to transform all 
the material information, v/hich is related to principle mate- 
rial axes, to each elemental coordinate system. To trans- 
form the material orientation given in terms of a global x - y 
coordinate system into each elemental coordinate, the follow- 
ing preliminary calculation was used. 

Given nodal positions of an element, x and y, then 
the transformation angle c() is 


(j) = Tan 


-1 


Y • “ y • ) 

13 _A 


X . - X. 

I 3 1 


as shown in Figure 1-17 , 


where, subscripts i, j correspond to the nodes 
numbered I and J, respectively. Since we know the 
global fiber direction angle a, 

3 = a - <() . 


These transformations were performed for 345 elements 
and 3 fiber orientations (45°/90°/-45° — the same as the 
photoelastic specimen, Ph - C type) by a preprocessor. This 
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X “ Y Global Coordinates 
Xe-Ye Elemental Coordinates 


Figure 1-17 

Transformation Diagram from Elemental to Global Coordinates 
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data preparation has now been completed and checked. 

The variable contact area between pin and hole, as load 
is changed, v/as simulated by using constraint equations at 
particular nodes on the loaded boundary of the hole. In our 
first computation/correlation attempt, eight nodes out of 
a total of nineteen around the half circle of the hole with 
which the pin could make contact Were constrained to have 
free displacement in the circumferential direction only. 

That corresponds to a contact angle of 80 degrees. All the 
nodes along the axis of symmetry were assumed to be simply 
supported. Along the loaded edge of the model, a uniform 
displacement field was used as a forced boundary condition. 
The value of the displacement imposed was 0.02 inch in the 
y direction (parallel to the axis of the symmetry) . 

The output from the main analysis program is also in 
terms of element coordinate systems. Such output includes 
three principal stresses (maximum and minimum normal stresses 
and maximum shearing stresses) at the centroid of the tri- 
angular elements and three stresses along the material axes 
(longitudinal, transverse and shearing) on each layer. 

The following procedure was used to check the consis- 
tency of the equations used in the program. 


1) Using the stress output along the material axes 
and the strains in terms of the same coordinates. 


VJher e , 
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and 


[Dol = 


X liT 


-X 0 
0 0 ^ 


LT 


again, where the subscripts L and T represent 
the fiber and transverse directions, respec- 
tively. But from the general strain tensor. 


[£ol = 




LT 




The characteristic equation for the eigen func- 
tions which determine principal strains is 


= 0 


from which the eigen values are 

X = 1 




2) Accordingly, the direction of the maximum value 
ip is 


ip = Tan 


- 1 


-^(e - £') 

Ylt ^ 


where, obviously, in this case, 

niax max 

and the maximiim shear strain y is 

'max 


'max max min 


3) Using the laminated plate theory and strains 
calculated above, principal stresses were com- 
puted and compared with those from the program 
output. The results were identical. 
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As noted earlier in this section of the report, strain 
contours are the principal means of correlating with the 
photoelastic results. To plot the strain contours, the fol- 
lowing procedure was used. . The maximum shearing strains at 
the elements' centroids were determined using the calcula- 
tion method described above. These centroidal values were 
then used to compute values at the nodes, by a weighted 
least-squares fit. The assumed fitting function was the 
second order, two dimensional polynomial 

P(x,y). = a^ + a^x + a^y + a^x^ + a^xy + a^y^ 
and the weighting function was. 


’w^(x,y) = 


(x,y) 


where f^(x,y) is the centroidal value of the element i, and 


r = /(x - x^) 2 + (y - y^) 2. 

The least square approximation technique to get the constants 
a^ (where i = 0; ••• 5) was checked with a one-dimensional case 
on the IMLAC graphics terminal. The resulting agreement was 
excellent; the curve passed through all the arbitrarily gen- 
erated data points and was smooth in the regions outside of 
the end points. Furthermore, no extra calculations such as 
an extrapolation were needed to calculate values at the nodes 
on the edge of the model. After getting the nodal values of 
strain, contour lines were drawn on the graphics terminal 
screen, by applying a linear interpolation method locally, 
based on the three nodal values of each element. 
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A first sample calculation was completed during this 
period, and the plot shown in Figure 1-18 exhibits patterns 
of maximum shearing strain distribution encouragingly similar 
to those of the photoelastic fringes,. The photoelastic pic- 
tures show half fringes (Figure I-19al and full fringes (Fig- 
ure I-19b). . The fringe number 1 is equivalent to a maximum 
shearing strain value 

Y - 0.18917 X 10 ^ inch/inch, 
max 

This is the same scaling factor used in the contour plotting 

of analysis results. Thus, a fringe number 2 corresponds to 

_ 2 

a value of 0.37834 x 10 inch/inch. The half fringe 

number 1 on the photoelastic picture in Figure I-19a is equiv 
alent to = 0.09458 x 10 ^ inch/inch. The load applied 

in the photoelastic experiments was 1500 pounds, and the 
post-calculated load of the finite element model was 1480 
pounds. 

Although the fringe patterns are similar, the finite 
element analysis results show almost exactly twice as many 
fringes as the photoelastic pictures show. The reasons for 
this are not yet clear. The somewhat different shape of the 
fringes may be caused by one or more of the following reasons 

1) Possible Material Differences. The material 
information used in the finite element analy- 
sis may differ from the characteristics of 
that used in photoelastic tests. Careful 
coupon testing would be required to resolve 
this; "sensitivity checks" using arbitrary 
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Figure 1-18 

Contour Plot of Maximum Shear Strain 
from the Finite Element Analysis 
(applied load, 1480 lbs.) 
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variations of material properties in the anal- 
ysis, however, may be worthwhile. 

2) Contact Angle. The contact angle correspond- 
ing to the analytical contour plot has not 
been iterated. It would be extremely fortui- 
tous if the assumed contact angle of 80° was 
correct. It is the purpose of this research 
to arrive at the correct contact angle by 
iteration, and this will be done. 

3) Slightly Different Model Size. The width of 
the analysis model is 2.4 inches and for the 
photoelastic specimen 2.48 inches. The dis- 
tance from the center of the hole to the free 
end is 2.0 inches for analysis model and 2.35 
inches for the specimen plate; thickness was 
the same for both models. 

4) Other Experimental Inaccuracies. Calibration 
of photoelastic coating, effect of grips, etc. 

D. Plans for Upcoming Period 

For the immediate fur ture, we will attempt to verify the 
material properties of the model and explore other potential 
sources of discrepancy between the fringes in the photo- 
elastic results and those of our finite element analysis. 
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CAPGLIDE (Composite Aircraft Program Glider) 

(R. J. Dienfendorf, H. J. Hagerup, G. Helwig) 

FAA certification in the "experimental" category was ob- 
tained for the first all-composite ultra^light glider, the 
RP-1, on August 25, 1980. A series of increasingly extended 
automobile tows from turf was then conducted at Schenectady 
County Airport, Schenectady, N. Y. during late July and 
August. The first public flight demonstration took place on 
September 16, 1980 with extensive media coverage. From first 
lift-off to public demonstration, the aircraft has made a 
total of sixteen successful flights using automobile-tow with 
release altitudes up to 700 feet. Significant gain in alti- 
tude above the release point was observed on at least one 
flight, although no attempt was made in these early trials 
to achieve soaring flight. Airplane traffic at the airport 
caused the tower to restrict all RP-1 flights to a relatively 
small area; this prevented a search for thermals. Mr. M. 
McCarron, a professional pilot, has piloted all flights to 
date. His reaction to the aircraft has been very positive, 
emphasizing its good inherent stability characteristics and re- 
sponsive control . Photographs of these activities at Schenec- 
tady County Airport are reproduced in Figures II-l through 6. 

The only untoward incident was experienced when the 
attempt to launch the fourteenth flight was abandoned, and a 
relatively mild, partial ground loop followed which produced 
minor structural damage to the skid assembly. We believe 









Figure II-6 
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that the difficulty was due to the use of an excessively 
elastic tow rope and too-rapid tow acceleration technique. 
Calculations confirmed that this caused the cable-aircraft 
combination to act somewhat like a slingshot, causing the 
glider to surge forward following the abrupt drop of friction 
from static to sliding during the first acceleration of the 
automobile. Then, after the potential energy in the cable 
was exhausted and tension in the cable dropped to zero, the 
glider momentarily slowed down. The pilot, sensing yaw and 
a reduction in directional control at that point, released 
from the tow rope. The aircraft yawed left, coming to rest 
in an area of sparse tall grass, which was probably the di- 
rect cause of the final ground loop. The damage to the skid 
assembly was quickly repaired, however, and the flight test 
program continued with a conscious effort to accelerate the 
tow car quickly only after the glider was seen to move. 

Operations are presently being moved to another site 
where winch tow to higher release altitudes will be possible, 
and there will be fewer restrictions on operating areas. 

The RP-1 landing gear is being modified to incorporate a 
wheel, thus facilitating launch by winch-tow. Release alti- 
tudes up to 1400 feet are expected. This second phase of the 
flight test program is designed to obtain performance data 
such as stall speed and glide polar for the RP-1. These tests 
can be performed safely only at considerable altitude. 

Meanwhile, work on the conceptual design and the pre- 
liminary design of the second aircraft has started. It now 
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appears that the RP-2 will be a 44-foot span tapered wing, 
ultra-light, all-composite sailplane with a fully-enclosed 
conventional fuselage. The aircraft, in order to ultimately 
have a self -launch capability, is being designed with space 
and structural provisions for an auxiliary engine. We ex- 
pect to conduct, however, initial test flights without an 
operational engine installed. The airfoil will be the thick 
version of the variable geometry airfoil developed by Boeing 
in cooperation with RPI and reported in the last semi-annual 
progress report. Preliminary experiments indicate that a 
structural design for the RP-2 wing, consisting of one-sided 
Kevlar-on-foam sandwich skins, will have the necessary 
strength and stiffness. This design approach allows a very 
simple fabrication technique. The one-sided sandwich skins 
can be laid up flat and cured and then bonded to the rib-spar 
shear-webs internal wing structure. The lower bending stiff- 
ness of the cured, one-sided sandwich skins allows them to 
take the shape of the ribs to which they are bonded. Wing- 
molds of the type used in the construction of the RP-1 will 
thus hot be necessary, saving time and effort. A preliminary 
three-view drawing of the RP-2 is shown in Figure II-7, and 
a preliminary version of the V-n diagram defining the loads 
envelope appears in Figure II-8. 

Activity in the Composites Fabrication Laboratory during 
the present reporting period has again been focused on com- 
posites structural techniques, with each student individually 
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or as a member of a team, undertaking a one-semester project. 
Generally, each project involves design, fabrication and 
testing. A partial listing of project topics selected for 
the Fall Semester follows: 

1 . Tapered graphite tubes 

2. Graphite-Kevlar hybrid leading edge 

3 . Variable camber wing section 

4. Fire-retarding composite materials 

5. Filament winding 

6. Shear test procedures 

7. Foreign object damage 

8 . Honeycomb sandwich 

9 . Aeroelastic box-beam wing 

Some of the projects are described in additional detail as 
follows; 

Project 1 will apply fabrication techniques and test 
data developed in the small-scale project work last year to 
construct a full-scale, circular cross-section, tapered fuse- 
lage boom. This is the type of construction which we antici- 
pate using for the aft portion of the RP-2 fuselage. The 
boom will consist of 0° and ±45° layers of graphite prepreg 
for bending and torsion stiffness respectively, the number of 
layers to be determined from deflection and stress calcula- 
tions for the boom when it is subjected to the design tail 
loads. The weight of the boom is expected to be approximately 
10 pounds. 

Project 2 explores the feasibility of replacing the 
present balsa wood leading edges by graphite-Kevlar hybrids. 
Kevlar cloth is laid up between two layers of 0° graphite 
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prepreg, and curing takes place in a mold incorporating the 
desired leading edge contour. The result should be a tough, 
light weight and easy- to- fabricate leading edge. 

Project 9 continues work from last year to predict 
coupled bending and twisting in a box-beam wing. As a re- 
sult of model construction and tests of such a wing last 
year, this coupling is now well understood and can be pre- 
dicted using the classical equations for orthotropic plates 
or finite element analysis. The next step is to apply the 
concept to a forward swept wing so as to increase bending- 
torsion divergence speed to an acceptable level without ex- 
tensive internal structure and/or added weight. 
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COMPAD (Computer Aided Design) 

(M. S. Shephard) 

The objective of the computer aided portion of the com- 
posites project is to provide computer "tools" for the analy- 
sis and design of composite structures. COMPAD ' s major 
thrust has been in the finite element area with effort di- 
rected toward implementing and improving the efficiency of 
existing analysis codes at RPI, developing our own inhouse 
analysis codes and developing interactive computer graphic 
preprocessing and postprocessing capabilities. The majority 
of recent effort has been in the development of the inter- 
active computer graphics programs because they represent the 
type of user interface that is required for effective use of 
advanced analysis techniques in the design of composites. 

That is, the importance of making new programs available is 
secondary to the effective use of those in operation at this 
time. 

The status of the RPI finite element software system is 
best indicated by the fact that several analysis programs 
(one inhouse program and several outside programs implemented 
on the RPI computer systems) , an interactive two-dimensional 
finite element preprocessor and an associated postprocessor 
exist and are being used more or less routinely. 

Several requests have been made to change the acronym of 
the software system from RPIFEP to something that can be ar- 
ticulated more easily. Therefore, the acronym POFES has been 
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adopted; this stands for People Oriented Finite Element 
Software. 

Effort during the current reporting period has been con- 
centrated in the following three areas: 

1. Implementation of a nodal resequencing algo- 
rithm for use in reducing finite element 
solution times. 

2. Development of an interactive graphic two- 
dimensional postprocessor. 

3. Development of an interactive graphic at- 
tribute editor for the insertion of material 
properties, applied loads and boundary condi- 
tions to the finite element models generated 
with the two-dimensional preprocessor. 

The first two areas above have been completed to the 
point that versions of the programs are available for general 
use, while the last item is in the early stages of develop- 
ment. Further details of the progress made in each of these 
areas are given in the following sections. 


1. Nodal Resequencinq Algorithm 

Efficient solution of the algebraic equations for a 
finite element model requires that the equation solver used 
accounts for the sparse and banded nature of the coefficient 
matrix. All general purpose finite element codes employ 
equation solvers that operate on the bandwidth^. 


Bandwidth = the maximum spread in rows or coliimns of non- 
zero elements on either side of the main diago- 
nal of the matrix . 


1 
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wavefront^ or skyline of the coefficient matrix. The 
solution time required for a model is more or less propor- 
tional to the square of the bandwidth, wavefront or skyline. 
The size of these parameters is dependent on the ordering of 
the unknowns in the coefficient matrix which is, in turn, 
defined by the ordering of the node or element niimbers of 
the finite element model. For most problems, excepting those 
with unusually regular geometries, the numbering of nodes or 
elements that yields maximvim solution efficiency is not ob- 
vious. Therefore, it is desirable to have an automatic 
means for ordering the node or element numbers that will 
minimize solution time. 

Several approaches have been taken in the development 
of reordering algorithms, with those based on graph theory 
being the most efficient and versatile. The algorithm chosen 
for implementation in POFES was developed by Gibbs, Poole and 
Stockmeyer^^^. The implementation of this algorithm is gen- 
eral, so it can be used to renumber node numbers or element 
nvimbers, depending on the finite element code chosen to ana- 
lyze the problem. This program runs very quickly and has 


^Wavefront = the profile of first non-zero elements, view- 
ing the matrix from left to right and row by 
row. 

^^Skyline = the profile of first non-zero elements, view- 
ing the matrix from top to bottom and column 
by column. 

^^^Gibbs, N. E., W. G. Poole, Jr. and P. K. Stockmeyer, "An 
Algorithm for Reducing Bandwidth and Profile of a Sparse 
Matrix", Siam J. on Niam. Analysis, Vol. 13, No. 2, April 
1977, pp. 236-250. 
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been found to result in large time savings in all but very 
small problems. 

2 ♦ Two-Dimensional Postprocessor 

The function of this postprocessor is to make it easier 
for the user to obtain a graphical representation of the re- 
sults of a finite element analysis. The basic capabilities 
of the postprocessor include deformed shapes, contour plots 
of various stress components and the ability to interactively 
obtain selected numerical output. In addition, the. user has 
the option to have the results of the problem sent to the 
line printer. 

The postprocessor has been designed to interface with 
the two-dimensional preprocessor to obtain the basic model 
definition and to interface with several analysis programs 
for their output. At the current time the postprocessor is 
only linked to our inhouse program FEPROQ, but it has been 
designed so that it can also be interfaced to other analysis 
programs such as SPAR, SAP-IV and ANSYS. 

The way we plan to use the postprocessor is best ex- 
plained by an example. Consider a bracket of uniform thick- 
ness (Figure III-l) fixed at the left end and loaded along 
the right-hand edges of both the hole and hook. Figure III-2 
shows a close-up of the right-hand portion of the bracket 
with the undeformed mesh in dashed lines and the displaced 
mesh (with a scale factor of six applied) in solid lines. 

In addition to displaying the displaced shape in terms of the 
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Deformed Mesh Superimposed on Undeformed Mesh 
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meshes, the deformed and undefomned outlines of the object 
can be shown. The user can also interactively change the 
scale factor applied to the deformed shape. 

The numerical values of the displacements and stresses 
at the node points of the element mesh are also available to 
the user interactively. Figure III-3 demonstrates the use 
of this option. This figure shows a close-up of the right- 
hand portion of the bracket with a deformed mesh superimposed 
on an undeformed outline. This plot indicates that the maxi- 
mum displacements occur at the bottom of the hook. Thus the 
user has interactively selected (by pointing with the light- 
pen which placed a "box" around it) the node that has the 
numerical values of displacement and stress which are de- 
sired. These values are printed to the far right under the 
associated node number for that point. 

To obtain more infomnation on the stress field through- 
out an object, the user can view contour plots of the maxi- 
mum principle stresses (Figure III-4) , minimum principle 
stresses (Figure III—5) and maximum shear stresses (Figure 
III-6) . The user can specify the contour level values in 
three ways; (1) the specific values of stress levels to be 
contoured, (2) a number of desired levels (Here the values 
for that number of equally spaced contours and their locations 
will be shown.) or (3) a starting value and increments of x 
and y coordinates and have the other levels determined at 
each point, incrementing in both directions till the maximum 
and minimum values of x and y on the body are passed. 
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Figure III-3 

Obtaining Selective Numerical Output 
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Figure III-6 
Maximum Shear Stresses 
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Although these contour plots indicate the basic trends, 
they can be enhanced to yield more information. Figure III-7 
demonstrates two features available for doing this. The 
first is the labels that have been placed on the contour 
levels. The values of the contours are shown in the far left 
portion of the figure. The second feature is the ability to 
highlight a contour level by pointing at it. This has been 
done for contour 11 in Figure III-7. 

Details of the means by which contours are obtained may 
be described as follows: 

a) Starting with the set of nodal values describ- 
ing an element, the element is broken up into 
a set of triangles. (The three-noded element 
is represented as a single triangle, the four- 
noded element is represented as four trian- 
gles, the six-noded element is represented by 
36 or 64 triangles and the eight-noded element 
is represented by 64 or 128 triangles.) 

b) The values of the stress component being con- 
toured are calculated at each of the trian- 
angles' vertices on the element. (The above 
two operations are efficiently accomplished 
in the natural coordinate systems used in the 
finite element representations of the ele- 
ments . ) 

c) Positions of portions of the desired contours 
of the chosen stress are then determined as 
intersection points of that value with the 
triangle's edges. The locations of these in- 
tersections are based on a linear interpola- 
tion of the corner values. Each of the 
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detected intersections is then properly con- 
nected to the intersection of that value on 
one of the other two edges. Since the three 
corner values of a triangle uniquely define 
a plane in terms of the value of stress being 
contoured, there is no problem determining 
the proper connections of intersections. 

One must exercise caution when determining the values 
of nodal stresses used by the contouring algorithm. It is 
well known that the direct evaluation of nodal stresses, by 
use of stress-displacement matrices, yields very poor stress 
values. In addition, the element by element evaluation of 
the nodal stresses will result in discontinuous contours at 
the element boundaries because the values at the nodes along 
an edge will be different for the two elements. The second 
problem can be overcome by a simple averaging of the nodal 
values for each element sharing a node. This averaging 
would yield continuous but inaccurate contours, due to the 
first problem. 

It is also well known that the most accurate values of 
the stresses are obtained by evaluating them at the locations 
of the niimerical integration points. However , the only 
points that prove convenient for use by a contouring algo- 
rithm are the node points. Therefore, a method for obtaining 
accurate nodal values based on the values at the numerical 
integration points is needed. There are two satisfactory 
means of doing this. The first and most accurate is referred 
to as "global stress smoothing" . The computations required 
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to do this, however, are more extensive than the original 

finite element analysis and cannot be justified for an inter- 

* 

active program. The second method is sufficiently accurate 
but requires far less computation. In this procedure the 
stresses are first smoothed on the elemental level using a 
least squares fit, based to a shape function one order less 
than that used for the displacement field. Averaging is 
then performed for nodal values using all the elements meet- 
ing at that node. 

3. Attribute Editor for The Two-Dimensional Preprocessor 

The preprocessor discussed in the last two progress re- 
ports gives the analyst the ability to define finite element 
representations for any two-dimensional geometry. However, 
before a finite element analysis can actually be performed , 
information about the material properties, loads and boundary 
conditions must be specified for the object. Since this in- 
formation is dependent on the class of problem, a seperate 
program (referred to as an "attribute editor") is required 
for each class. Effort is currently underway to develope an 
interactive attribute editor for stress analysis problems. 

The operation of the attribute editor is controlled by 
a main control page (Figure III-8) that indicates the basic 


Hinton, E. and J. S. Campbell, "Local and Global Smoothing 
of Discontinuous Finite Element Functions Using p. Least 
Squares Method", Int. J. Num. Meth. Engng., Vol. 8, 1974, 
pp. 461-480. 
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Figure III-8 

Attribute Editor Control Page 
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program modules. Thus far the modules that have been started 
include the point load, boundary condition and material prop- 
erty editors. 

The function of the material property editor is to give 
the user a means for specifying the material constants asso- 
ciated with the elements of the various subregions. Figure 
III-9 shows the use of the material property editor in the 
application of properties to a deep beam made-up of two dif- 
ferent isotropic materials. This editor is being designed 
to allow composite materials' properties to be defined con- 
veniently. 

Boundary conditions are specified by means of the menu 
page shown in Figure III-IO. The user can use this menu to 
interactively fix any of the possible translations or rota- 
tions (in bending problems) of the nodes; fixing one node at 
a time along a portion of an edge curve or in an entire mesh 
subregion. For the deep-beam example, the user has fixed 
both translations at the bottom-left corner and the vertical 
motion of the node next to that. He has also fixed the two 
nodes on the right edge of the bottom to simulate the beam 
sitting on bearing pads at either end. 

The point load editor operates in much the same manner 
as the boundary condition editor with the only difference 
being the specification of the load value. Figure III-ll 
shows the display given after the application of point loads 
along a portion of the top of the beam. 
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Figure III-9 

Application of Material Properties 
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Figure III-IO 

Application of Boundary Conditions 
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Plans for the upcoitiming period call for the major effort 
to be placed on the stress analysis attribute editor. The 
portions requiring the most work are the material property 
editor, edge-load editor and surface-load editor. In addi- 
tion, a program will be written to link the attribute editor 
to FEPROQ. 

Although the postprocessor in. its present form provides 
a useful capability, some additional features are desirable. 
One of these features which, however, will require substan- 
tial effort is the ability to properly contour stresses in 
the presence of material interfaces. Work is expected to be 
done in this area as well as on the attribute editor. 
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PART IV 

INSURE (Innovative and Supporting Research) 


A. Advanced Structural Analysis Methods for Composites , E. 
J. Brunei le. 

B. Ultrasonic Non-Destructive Testing of Composite Struc- 
tures , P. Das and H. F. Tiersten. 

C. High Modulus Graphite Fiber Surface Modification for Im- 
proved Interaction with Matrices , R. J. Diefendorf. 

D. Fatigue in Composite Materials , E. Krempl. 

E. Acoustic Emission of Composite Materials Undergoing 
Fatigue , H. Otsuka and H. A. Scarton. 

F. Resin Matrix Dependent Characteristics of Composites , 

S. S. Sternstein. 




IV-A ADVANCED STRUCTURAL ANALYSIS METHODS FOR COMPOSITES 


Senior Investigator: E. J. Brunelle 

The filamentary composite lamina is the basic building 
block of almost all composite structures . Metaphorically it 
contains the genetic code that governs the general behavior 
of laminate assemblages. The specialized behavior of a par- 
ticular composite structure Cits deformation characteristics, 
etc.) is determined by its cellular structure (the stacking 
sequence) and by its phylum (beam, plate, shell, etc.). It 
is suggested that to completely understand the behavior of 
various composite structures one must back-track through the 
path of its creation and first understand the behavior of the 
individual lamina in as illuminating and as efficient a repre 
sentation as possible. The correct representation, it is 
held, comes from two viev-?points, not one. The first view- 
point is that one must clearly understand the behavior of the 
lamina's constitutive constants under coordinate rotation 
since reflection and translation yield little useful informa- 
tion. Once this part of the representation is fully under- 
stood then one may move to the second viewpoint which holds 
that if the changes in response of the lamina to the syste- 
matic changes of external stimuli are understood (say by a 
very efficient and concise method of generating solutions to 
various stimuli ) then one may claim to possess a valid repre 
sentation of the behavior of a single lamina. To these ends, 
two broad avenues of investigation are being pursued. 
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1. Status 


The major thrusts of the last two reporting periods (pre- 
vious periods were not devoted to INSURE) were (a) the deline- 
ation of the rotational elastic properties of an orthotropic 
lamina, the key parameter being the nondimensional number 
B = which predicts "regular" or "anomalous" rotational 

properties depending on whether B > 1, or B < 1 respectively. 
Extensions of this work to composite structures and its ex- 
pected important role in off^design structural conditions was 
outlined and temporarily set aside to concentrate (b) on find- 
ing the most efficient representation of thin orthotropic 
structures and to find concise solutions for said structures. 
This work led to the useful concept of multiple affine trans- 
formations, several nondimensional numbers (the generalized 


D 2 + 2D g 

rigidity D* = — — — being the most prominent one for 

’^° 11°22 

plate structural solutions) and the existence of solution 
functional dependences for many classes of composite struc- 
tures problems. The usual form of these dependences was 




a n 
o 


b m 
o 


were a^/b^ is a length ratio in the affinely stretch- 


ed coordinate system, m and n are positive integers, and p is 
a plus or minus, integer. Finally (c) , a systematic effort 
was undertaken (and is still in progress) to provide generic 
solutions for a wide variety of problems . The largest body 
of solutions to date presents elastic buckling coefficients 
of thin plates for a wide variety of loadings and boundaries; 
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the only parameter in most of these solutions being D*, which 
takes values in the range o < D* < 1. 


2. Progress During Report Period 

Two more closed form exact solutions have been obtained. 
The first of these solutions yields buckling coefficients and 
vibration frequencies of uniaxially loaded simply supported 
rectangular plates bonded to a core material. The core is 
modelled as a Winkler Foundation (its inertial effects are not 
accounted for) thus the frequency expressions are valid only 
for the low-to-middle range of the usual classical vibration 
domain; for most practical purposes this will not be a prob- 
lem. Both solutions come from the homogeneous partial differ- 


ential equation (in affinely stretched y^ and t^ coordi- 


nates) ; 
3 ‘‘w 


3x 

o 


+ 2D*- 


3^w 


3x ^dy ^ 
o -^o 




3 ^w 


+ kw + 


v/D 


11 '^o' 


3 ^w 
3t_^ 


= 0 


where 


3t 


= ph- 


3t' 


and k is the foundation (core) modulus. 


Defining a buckling coefficient k^ = 


'^° 11°22 


and letting 


w(x^, Yq) assvuae the usual multiple half sine wave behavior in 
both the x^ and y^ direction (the time variable is suppressed) 
we find that 


k - 2D* = 
o 


0 

g 

V . 

2 

4. 

f 

a 

o 

2 

1 + ^ 


4 

r°J 


mb 

o 

^ J 


°22 




and n = 1 with the 


following detailed information. The minimum k^ - 


2D* values 
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k fhY'* 

are given by (k - 2D*),,. = 2/^ where K* = 1 + •=; — - — 

^ o Min. D 22 l^J 

and the corresponding values of a^/b^ are given by 


> 1 , 


^^o'^^o^Min. (k - 2D*) 
o 


where m is the mmber of half- 


waves in the X direction. Furthermore, the values of a /b 
o o o 

at which the buckling coefficients are equal are given by 


j. 1 

O' o m = m + 1 


'm(m + 1] 


-. These results are presented 


in Figures IV-A-1, -2 and -3. The dominant role of the core 
as it becomes progressively stiffer is quite evident. When 
the time variable is not suppressed and we let 

w(x , y , t ) = A sin — — sin ° e ° a solution for the 

0-0 o a D 

o o 


non-dimensional frequency Q, 


mno hr D 


) « 'i 2 


is found by utiliz- 


ing the construction k^ = 2 (D* + ^ where o £? £ 1 (if only 

compressive loads are assumed). The result is 


+ 1 - K* - n'^ /n'^K* = — j 


K* na 


** r* r * v.lf mb 

+ 2 ^ . -2^ + 1 i 

K*'^na 

V. y 


D* E 

Noting that the maximum values of and are unxty, the 

n^ 

coefficient they describe is very limited, thus the left hand 

mb 

side may be plotted versus — — — — with the limiting combina- 

K*'*na 

o 

tions of and — = (0, 0), (0, 1), (1, 0), (1, 1). 

i/K* n^ 

Some typical results obtained using the developments dis- 
cussed above are shown in Figures IV-A-4 to -9. It should be re- 


iterated that these graphs represent generic solutions; that 
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Figure IV-A-6 


Stability Coefficient vs. As 
pect Ratio of an Orthotropic 
Plate Clamped on All Edges an 
Compressed Between One Pair 
Opposite Edges 










Figure IV-A-9 

Frequency Parameter vs. Stability Co 
efficient of a Circular Orthotropic 
Plate with Hinged Edge. Under Com- 
pression „ - .2 
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is, they are valid for any material and any dimension ratios 
for the specified rectangular or circular geometry in affine- 
ly stretched coordinates. During the process of checking 
these results with the existing particular solutions in the 
literature it was realized that by proper manipulation these 
existing particular solutions can be converted into generic 
solutions, thus greatly extending the usefulness of existing 

solutions (and eliminating the need to duplicate the corres- 

( 1 )* 

ponding generic results). For example Hearmon presents an 
energy solution for the center deflection u of a clamped 
plate under a uniform distributed load q as; 


q - 

292u 


D 


11 


D 


+ 


22 , 4 ^^^12 + 
7 


^ a"* b"* ' a^b" 

A little algebra produces the generic result 


292u 7 


f *N 

a 

2 

fb 1 

o 


o 

b 


a 

o 

V. J 


oJ 


V2 



which is just a single curve in the variable a^/b^. Note now 
much more useful the second form is than the first form. 

The second exact solution obtained in this period yields 
moisture diffusion results for a specially orthotropic slab. 
These results have been derived for the x^ = o,c faces at 
different equilibrium moisture concentrations at t = o and 
at t oo; i.e. , the upper face goes from a moisture concen- 
tration r° at time zero to r” at steady state (t ^ <») . In a 
u u 


Bracketed numbers in this text refer to references on page 
115 
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like fashion, the lower face goes from a moisture concentra- 

00 

tion r° at time zero to at steady state (t <=°) . This 
allows the more typical condition of gases of different mois- 
ture concentrations to exist on each side of the slab' and 
yields the essential case of a liquid/moist gas condition for 
treating submersible structures. It is easily shown that the 
one-dimensional and the three-dimensional solutions are re- 
lated by multiplicative factors; long-time and short-time 
solutions are easily constructed, and an integration over the 
slab volume yields the total moisture content G(t^) as a func- 
tion of non-dimensional time t^. As is well known, the short- 
time (three-dimensional) solution for G(t) is linear in the 
square root of elapsed time [see, for example (2)]. In the 

variable t one has the counterpart relation G(t ) = m/t~ for 
o o o 

4 

short-times. The slope m = — 


c 

a 


1 + 

1 + ^ 


a 

b 


o 

o 

V J 






is linear in 


both the thickness ratio — and the aspect ratio a /b , where 

a " o o 

o 

a , b and c are the slab dimensions in the affine space 
o o o ^ 

X , y , z and t = t/c^^ (see Figure IV-A-10) . The affine 
o o o o o ^ 

transformations used in the moisture diffusion problem being 
Xi = where are the moisture diffusion coefficients 

and the slab dimensions in real and affinely transformed 
space are related by: 


a^ = a/v/K^ , 


bo = b/^ 


and c = c//kT. 
o ' 3 


Some brief solution details are given. Firstly, the equation 


governing moisture concentration F is given by K. 


3^r 


•9Xi^ 


II 

3t 
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which is affinely transformed into Xj_ space to yield, 

_ ar 

9 Xi 9 Xi (t will be transformed to a suitable t^ later) . 

A solution is formed from the steady state solution as t 

T 

and a transient solution. Thus, T (Xj^/t) = ^ . 

00 CO ^3 

From previous comments it is realized that = (T^ - r^) — + 

^o 

oo T 

r^, and one forms a separable representation of F as 

(xj^) T2 (X2) ^3 (X3) T (t) , where the value of F^ on all six 

T 

faces must be zero at all times, and as t -> <=°, F o through- 

T 

out the entire space. When t = o, F^Cx^) “ ^00^X3) + i (Xj^/O)/ 
o 0^3 o 

where F = (F - F„)^ + F.. Insuring these requirements 

O U 36 C 36 

O 

leads to the desired result. 


00 00 00 3 

r(Xi, X2^ X3vtd = + (r^ " ^ 

O 


q - ^2 ^3 

■ ■ ■■sin (2m - 1) ir — sin (2n - 1) ir— sin stt— . 
m=l n=l s=l(2m - 1) (2n - 1) (2s) a b c 

o o o 


00 00 00 

I, X, X. 


.exp(-l) 


(2m - 1)^ + (2n - 1) ^ 






f ^ 

a 

2 


0 


+ s2 

b 



oj 

— 



TT^t 


where t = t/c f = F° - F°° and g = F° - F?- 
o o u u ^ £ I 

An integration over the slab volume yields the total moisture 

content G(t ) = /F (x • / t )dV, from whence the short time solu- 
o Y ^1 o 

tion can be deduced. Also a one-dimensional solution for G(t ) 

o 

can easily be deduced from the three-dimensional result. It is 


G(t ) - JsO’° + r°) 

o u £ 

t oo CO 0 0 — 

<''u " ''F - 


e - (2s - 1) ^TT^t 
8 “ o 

tt 2 (2s - 1)2 


1 
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3. Plans for the Upcoming Period 


Two items appear to be logical areas for investigation. 

The first is to continue solving problems of importance which 
possess specially orthotropic constitutive relations . The 
second is the difficult one of trying to find the most effi- 
cient representation of solutions for problems with generally 
orthotropic constitutive relations. One remembers that a 
specially orthotropic plate has only four elastic constants 
whereas the generally orthotropic plate has six elastic con- 
stants and has its behavior complicated by the twist-coupling 
associated with and D 2 g. The best generalization that one 

would dare hope for is that only one new "cluster constant" 
will appear, thus making generic solutions depend only on 
this new "cluster constant" and D*. From past R.P.I. research, 
however, one has an excellent idea of how solutions vary with 
D*. Therefore, it should be possible to concentrate on the 
solution variations with respect to just the (one) new cluster 
constant. If this speculation proves to be valid, then one 
should be able to predict directly the generic effects of 
changing individual plies upon the solution of a symmetric 
composite plate. This would appear to be a significant step 
forward in understanding the behavior of a large class of com- 
posite structures. In the upcoming period, this objective 
will be emphasized. 
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IV-B ULTRASONIC NON-DESTRUCTIVE TESTING OF COMPOSITE STRUCTURES 

Senior Investigators: H. F. Tiersten 

P . Das 

Ultrasonic waves are being expolited both in research and, 
increasingly, in routine non-destructive testing as a prom- 
ising means of discovering and assessing manufacturing flaws 
and load-induced damage in composites. Our efforts in ultra- 
sonic NDE are three-fold: improving image quality through 

data processing and display techniques, transducer development 
and material flaw image identification. 

1. Status 

Improvements in ultrasonic image quality have been ob- 
tained through use of a) a high resolution display monitor 
and b) color -coding of the intensity data. Both experimen- 
tal and theoretical progress has been made in improving 
trapped energy mode mosaic ultrasonic transducers, including 
Fresnel lens types. Support of other project activity, using 
ultrasonic NDE is being provided on a service basis. Examples 
to date include delamination patterns around holes in bearing 
test specimens. 

2. Progress During Report Period 

To obtain high quality color-coded images, the raw stored 
data is now being routinely transferred to the Image Pro^ 
cessing Laboratory which has a DEANZA image processor with 
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a PRI^^E host computer. The transferred data is being stored 
in the PRIME file system. For display or processing, the 
data is then transferred to the image refresh memory. The 
color coding of the data can now be achieved through a pro- 
grammable intensity transformation table. The image is dis- 
played in color on a high resolution monitor with zoom and 
scroll capability. 

A series of ultrasound scans have been performed as a 
"blind test" , i. e. , a test in which the expected results of 
the scans were not known to the experimenter. The graphite 
composite samples were prepared and supplied by Professor S. 
Sternstein of the Materials Department at R.P.I., and this 
report presents the results which were obtained. 

Six samples were used in the test. Each sample was a 
circular disk about 4 . 3 cm in diameter and .18 cm thick, 
with a mass of 2.5 gm (Figure IV-B-1) . Since the samples 
appeared to be identical to the unaided eye, they were each 
identified by a code written on with a grease pencil. The 
codes were lA, IB, 2A, 2B, 3A and 3B. Also, an arrow was 
drawn on each disk to show the proper orientation for the 
scan. Although the meaning of these codes was not revealed 
until after the test results had been presented, they will be 
explained here for clarity. 

The letter part of the code describes the physical con- 
struction of the sample while the number describes the degree 
of damage caused by stressing the sample. This is shown in 
Table IV-B-1. The samples were stressed prior to ultrasonic 
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inspection by supporting them around their circumference with 
a circular jig and applying a force directed perpendicular to 
the surface at the center of the disk. For ultrasonic test- 
ing, each of the samples was mounted in a plastic clamp that 
held the edge of the disk without stressing the sample. Fig- 
ures IV-B-2 to 7 show the displays of the images resulting 
from ultrasonic scanning and image processing. 

The degree of damage to each sample is readily apparent 
from the displays. Samples lA, IB, 2A and 2B show no damage. 
Samples 3A and 3B very clearly show a well-defined region of 
damage in the center of the disk, with the shape of the dam- 
aged area dependent on the construction of the disk. Also, 

3A shows alternating bands of green and red, indicating some 
form of damage along the fiber lines. The failure mechanism 
that caused this has not yet been determined . In addition, 
the sensitivity of the ultrasound system is such that on many 
of the displays it is possible to see the code that was 
written on the sample. 

The color coding scheme used shows heavily damaged areas 
in blue, while areas of lesser or no damage are shown in red 
and green, respectively. The coloring was chosen by examin- 
ing a computer generated histogram of the image. The histo- 
gram shows the distribution of intensities in the image, and 
the coloring was adjusted to bring out as much detail as pos- 


sible. 
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3. Plans for the Upcoming Period 


In the next period we anticipate further development of 
our signal processing techniques, perhaps to include deconvo- 
lution. Transducer advances will include application of 
higher frequency ultrasound, from 5 to 15 MHz, for better 
resolution and further development of phased array transdu- 
cers. 


4. Current Publications or Presentations by 
Professors Das and/or Tiersten on this Subject 


1. Das, P., S. Talley, R. Kraft, H. F. Tiersten and J. F. 
McDonald, "Ultrasonic Imaging Using Trapped Energy Mode 
Fresnel Lens Transducers, Acoustical Imaging, Vol. 9, 
Edited by Keith Y. Wang, Plenum Publishing Corp., 227 
West 17th Street, New York, N. Y. 10011, 1980. 

2. Webster, R. T. , P. Das and R. Werner, "Ultrasonic Imaging 
for Nondestructive Evaluation of Composite Material with 
Digital Image Enhancement", to be presented at the 1980 
Ultrasonic Symposium, Boston, MA, November 5-7, 1980. 


TABLE IV-B-1 

CODE FOR IDENTIFICATION OF "BLIND TEST SAMPLES" 


Code 


Construction 


A 0-90 degree fiber lay-up 

B 0-45 degree fiber lay-up 


Code 

1 

2 

3 


Damage 

no damage, unstressed 
stressed, but not to failure 
stressed to first failure 
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IV-C HIGH MODULUS GRAPHITE FIBER SURFACE MODIFICATION FOR 
IMPROVED INTERACTION WITH MATRICES 

Senior Investigator: R. J. Diefendorf 

The mechanical properties of composite materials are 
strongly influenced by the strength and behavior of the inter - 
facial region between the matrix and the fiber. Not only does 
the interfacial region determine whether an initiated crack 
will cause delamination or fiber breakage, it also will de- 
termine if a broken fiber will continue to reinforce. 

1 . Status 

( 1 ) * 

Goan and coworkers measured the short beam shear 
strength (SBSS) of composites using various high modulus gra- 
phite fibers and showed that the SBSS decreases with increas- 
ing fiber modulus. Moreover, for composites containing un- 
surface-treated fibers with moduli greater than 400 GPa, the 
SBSS drops to an unacceptably lov7 value of about 15 LIPa. 

This very low value of SBSS for a high modulus composite 
limits the application of these fibers in components where 
high shear loading may be present. While it can be argued 
that use of high modulus fibers may lead to composites with 
inherently lower SBSS, the high value of 125 MPa observed in 
boron fiber composites is direct evidence to the contrary. 


'k 

Numbers in parentheses in this section refer to the refer- 
ences listed on pages 154 and 155 
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( 2 ) 

Herrick observed that etching the surface of high 

modulus carbon fibers resulted in composites with a remarkable 

increase in SBSS. However, the trend of decreasing SBSS with 

increasing fiber modulus remained. Although the absolute 

magnitude of SBSS increased, unfortunately the tensile 

strength often decreased by as much as thirty percent. In 

( 3 ) 

Herrick's other related work, he concluded that surface 
functionality of the fiber was more important than the sur- 
face area in determining the SBSS of the resulting composite. 


More recent v/ork on rayon precursor carbon fibers by Scola 
(4) 

and Brooks is in direct conflict with those of Herrick et 


The purpose of this study is to elucidate the molecular 
nature and surface energetics of graphite fiber surfaces. 
This will be done so that the phenomena occurring at the 
matrix-fiber interfacial region may be characterized, and 
knowledge of these phenomena may be applied to improve com- 
posite materials. 


2. Progress During Report Period 

In this report the essential parameter used to study the 
surface energetics of graphite fibers was the wetting by vari- 
ous liquids. Wetting was evaluated by measuring the contact 
angle of various liquids in contact with the fibers. The im- 
mediate difficulty v/ith this method is inherent with the 
fiber surface; that is 
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a) it is difficult to measure with reasonable 
accuracy the diameter [6 - 14vi] or perimeter 
of fibers, considering the fact that the 
cross-sectional area is not always circular. 

Also the cylindrical surface area per unit 
length varies from fiber to fiber. 

b) The roughness and smoothness and polarity 
of fibers is strongly dependent on the con- 
ditions used in processing the fibers (in- 
cluding any final surface treatment ren- 
dered to the fiber surfaces) . 

It is well documented that the surface of high modulus 
fibers HMU, HMS carbon fibers and others with relatively low 
shear strength, have far smoother surfaces than low modulus 
(200 GPa) fibers with reasonably high shear strength (Figures 
IV-C-1 and -2) . 



Figure IV-C-1 

Scanning Electron Micros- 
copy (SEM) of Longitudi- 
nal Surface of Low Modu- 
lus Fiber (E = 200 GPa) . 
Magnification 2 0,00 Ox. 



Figure IV-C-2 

Transmission Electron 
Microscopy (TEM) of 
Longitudinal Surface of 
High Modulus Dogbone 
Fiber (E = 680 GPa) . 
Magnification 6 6, 70 Ox. 
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In order to overcome the above difficulties, highly ori- 
ented pyrolytic graphite (HOPG) - the same as graphite mono- 
chromators used for x-ray and neutron analysis - was used to 
simulate the surface of high modulus fibers. While the aver - 
age crystallite preferred orientation and perfection of high 
modulus carbon fibers is significantly lower than those of 
HOPG, TEM studies indicate that the outer-most layer of high 
modulus carbon fibers does not appear to be drastically dif- 
ferent from HOPG. A comparison of some pertinent properties 
for HOPG versus Thornel P-75S is given in Table IV-C-1. 

A subgoal of this study was to determine if HOPG could 
be used for simulation. The surface of HOPG (Figure IV-C-3) 
can be characterized as consisting of two major parts: a) 

the basal planes (0002) with low bonding ability, hence poor 
wetting and b) the very active edges with high surface ener- 
gies with comparatively excellent wetting ability (schemat- 
ically illustrated in Figure IV-C-4) . The disparity in the 
behavior of the two areas is shown vividly in the sequential 
pictures of a drop of distilled water on a relatively pure 
basal region (Figure IV-C-5) and on the basal area containing 
some edges (Figures IV-C-6a through -6d). 

The surface energetic of the basal planes of HOPG would 
be expected to be very similar to that of the surface of un- 
treated high modulus fibers, because of the similarity of the 
basal plane orientation of HOPG to that of the outermost 
layers of high modulus fiber. Figures IV-C-3 and -7. 



TABLE IV-C-1 

SOME PROPERTIES OF HOPG AND HIGH MODULUS GRAPHITE FIBERS 


Temp . Range 
Over Which 


Properties 

Units 

HOPG' 

Prop. Was 
Measured 

C 

o 

Filament 

diameter 

0 

A 

6.710 

Room 




Density 

g/cml 

2.266 

Room 

Thermal con- 
ductivity 

o 

W/mK 

1700. (a)** *** 

8. (c) 

2-340°C 

90-340°C 

Electrical re- 
sistivity 

-4 

ohm/cm 10 

0.4 (a) 

0.15 (c) 

1-300 

1-300 

Thermal expan- 
sion coefficient 

H 

1 

0 

Slightly 
negative (a) 
28 X 10"® (c) 

300-1100 

Tensile modulus 

GPa 

Oil =1060 

Room 

(17) 

Carbon assay 

% Carbon 

99.9996 



Temp. Range 
, ^ Over Which 

Thornel P-75S Prop. Was 

(Fiber Axis) Measured 


lOy 

2.00 

150 at 21°C 

5 : _____ 

-1.1 X lo”^ at 2r°C 

500 

99 _____ 


Bracketed numbers in this table refer to references on page 155 

ick 

Along the (a) axis of hexagonal lattice 

** * . 

Along the (c) axis of hexagonal lattice 
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Figure IV-C-7 

Model for High Modulus Graphite Fiber 
(after R. J. Diefendorf and E. Tokarsky) 
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a) Experimental Technique 

The surface of HOPG was carefully cleaved under a stereo- 
microscope to obtain a sufficiently large basal region on 
which the contact angle of resins, hardeners and various liq- 
uids of known and unknown polar and dispersive components 
could be measured by using a goniometer. The same experi- 
ments were replicated on the surface of polished polytetra- 
f luoroethylene (PTFE) . The polar and dispersive components 
of HOPG and PTFE were evaluated by determining the slope and 
intercepts of the Young -Fowkes relationship: 



where: Yt refers to the equilibrium liquid-air surface energy 

J-J 

of the contacting liquid, 0 is the measured contact angle and 

Y? and y^ refer to the polar and dispersive components of the 

surface energy of the experimental liquid. Typical values 

for the liquids appear in Table IV-C-2. Similarly, Yg 
(i 

Yg refer to the polar and the dispersive contributions of the 

o ^ 

solid surface energy Yg- For the solid Tg = Yg Yg and for the 
TD d 

liquids Yt “ Yt + Yt- a detailed derivation of this equa- 

J-i i-i Li 

(7) 

tion, see Kaeble . 

Finally, the critical surface energy, Y^ / of HOPG and 

(5 7) 

PTFE was evaluated using the Zisman plot ' , and by using 

linear regression and measured contact angles, the surface 
energy of various hardeners and resins was estimated. A 
typical Zisman plot is shown in Figure IV-C-8. 
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TABLE IV-C-2 


SURFACE, ENERGY OF 

REFERENCE 

LIQUIDS AT 

20°-25'’C 


Reference 

Liquids 

Contact 
Angle 
6 on PTFE 

Contact 
Angle 
0 on HOPG 

Y ^ 
’D 2 

cm 

cm 

Tt 2 

cm 

Water 

112 

60 

51 

21.75 

72.75 

Glycerol 

95 

50 

37 

26.4 

63.4 

Methylene 

iodide 

75 

17.5 

48.4 

2.4 

50.4 

Glycol 

83 

30 

28.9 

18.8 

47.7 


Yjj = Dispersive surface energy component 
Yp = Polar surface energy component 

Y^ = Total, dispersive surface energy component and polar 
surface energy component 


3 




Slope = “0.0252 

Xc = '9-® 


Coeff = 0.9982 
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b) Results 

By using the Young -Fowkes relationship again; 


Cos j9 + 1 = 



+ 



where stands for the reference liquid total, d and p are 
the dispersive and polar components, respectively, and L and 
S are the liquid and solid, respectively. By rearranging the 
equation slightly; 


Yj^(Cos6+ 1) 




A plot of 


with 


Yj^(Cos0+ 1) 


versus 


should be a straight line 


f d) ^ 


Yg and 

yP 


as intercept and slope, respectively. 


Such relations are shown in Figure IV-C-9. The data shown in 
Figure IV-C-10 was used to compile the information contained 
in Tables IV-C-3 and -4 . A compendium of such data for a 
very wide range of resins, hardeners and liquids of interest 
in composites fabrication is presented in Table IV-C-5. These 
values were obtained with the methods reported here and from 
the literature. A summary of the surface energy of the solids 
PTFE and HOPG is given in Table IV-C-6 and compared with 
values in the literature. 



(l + Cose)/2(vd )'/2 
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Figure IV-C-9 

Determination of the Surface Energy Component of 
HOPG and PTFE at 20°C 



Cos Q 



7l 

Figure IV-C-10 

Zisman Wetability Plot of Various Liquids 
^on I-IOPG at 20 - 25°C 
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TABLE IV-C-3 

RESULTING DATA FOR SURFACE ENERGY COMPONENTS 
ON HOPG AND PTFE 


Reference 

Liquids 



PTFE 


HOPG 


Yt ( 1 + Cos 9) 

ij 


Y^(l + Cos 6) 


d 

h 


^L 

V. J 

2 

^L 

s J 


Water 

2.3448 

11.7673 

4.8778 

Glycerol 

0.8447 

8.5613 

4.7572 

Glycol 

0.8066 

8.2756 

____ 

Methylene 

iodide 

0.2227 

7.1330 

4.5959 


(Y 


(y 


PTFE' 


Y 


PTFE 


PTFE' 


Y 


PTFE 


Corr . Coef f . 


4.6034 

21.2 

0.1232 

0.015 

0.9499 

0.95 


= 21.2 ± 0.35 


(Y 


(Y 


HOPG 
d 


Y 


HOPG 


HOPG 


Y 


HOPG 


Corr. COeff, 


= 6.6244 

= 43.9 
= 2.1908 

= 4.78 

= 0.9991 


T 

YpTFE 


T 

"^HOPG 


48.7 ± 3.2 
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TABLE IV-C^4 

ZISMAN EQUATION RESULTS FOR VARIOUS LIQUIDS ON PTFE AND HOPG 


PTFE HOPG 


Liquids 

Cos 0® 

cm 

Cos 0® 

Y erg 

^LV 2 
cm^ 

Water 

-0.3746 

72.75 

0.5 

72.75 

Glycerol 

-0.0872 

63.4 

0.6428 

63.4 

Acetic anhydride 

0.6428 

32.7 

0.9903 

32.7 

n-Butyric acid 

0.8192 

26.8 

— 

— 

Acetonitrile 

0.7547 

29.3 


— 

Aniline 




0.9063 

42.9 

Formic acid (9.76%) 

— — 



0.6157 

63.6 


Linear Regression of 

Cos 0 vs . 

^LV^ 


PTFE 

HOPG 

(erg/cm^) 

19 . 6 

33.5 

Slope 

-0.0252 

-0.0125 

Coef f . Corr . 

0.9982 

0.996 

Ay (erg/cm^) 

— 

20 
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TABLE IV-C-5 

CONTACT ANGLE AND SURFACE ENERGY ESTIMATE FOR VARIOUS 
RESINS, HARDENERS AND LIQUIDS 




Contact Angle 
(Degrees) 

HGPG PTFE 

erg 

2 

cm 

Estimate 
for PTFE 

Y erg 
2 

cm 

PTFE 

Corrected 

V 

''l 2 

cm 

Estimate 
for HOPG 


RESINS 






1. 

Epon 828 

15 

85 

54.8 

49.1 

36 

2. 

Epon 871 

flexibilizing resin 

10-12 

67 

42 

37.6 

34.9 

3. 

Epon 152 
phenol Novolac 

25 

— 

— . — . 

— 

41 

4. 

Ciba MY-720 
Cresole Novolac 

32 

92 

59.6 

53.4 

45.7 

5, 

Ciba 179 

cycloaliphatic resin 

10-12 

78 

50 

44.8 

34.9 

6. 

Ciba 508 

<5 

70 

44 

39 

£33.8 

7. 

Tor Ion 400 (Amoco) 

20% imide in n-methyl 
pyrrol idone 

50 

99 

64 

57.4 

62.2 

8. 

Kelpoxy G-272-100 
epoxy terminated 
elastomeric copolymer 

30-31 

92 

59.6 

53,4 

44 


HARDENERS 






9. 

Nadic methyl anhydride 
(NMA) 

9 

50 

32.9 

29.5 

34,5 

10. 

Dodecenyl succinic 
anhydride ( DDSA) 

8 

50 

32.9 

29.5 

34.2 


Amide 






11. 

V-40 (Epon) 

22 

58 

37.3 

33.4 

39.3 


Amines 






12. 

Diethylene triamine 
(DTA) 

30 

80 

51,4 

46.1 

44.2 

13. 

Tetraethylene pent- 
amine (TPA) 

30 

80 

51.4 

46.1 

44.2 
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TABLE IV-C-5 (Continued) 


Contact Angle y. 

(Degrees) cm cm cm 

Estimate PTFE Estimate 
HOPG PTFE for PTFE Corrected for HOPG 


HARDENERS 
Coupling Agents 


A) Coordinate Types 

(titanate complex) 


14. 

KR44 (contains amine) 

31-32 

73 

46.7 

41.8 

— 

15. 

KR46^ 

spreads 

38 

27.17 

24.4 

<33.5 

16. 

KR55 

spreads 

51 

33.4 

29.9 

<33.5 


B) Chelate Types 






17. 

KR-112S 

8 

— 



34 

18. 

KR-134S 

18 

~ 

— — 


37.4 

19. 

Soap solution 

splash 

spreads 

— 



■ 

«33.5 


Oxidents 






20. 

Hydrogen peroxide 

30 


— 

— 

44 

21. 

Nitric acid (70%) 

15 




36.2 

22. 

Titaniimi trichloride 







(20% solution) 

55 

— 

— 

— 




Dispersives 






23. 

Acetic anhydride 

8 

50 

32.9 

29.5 

34.2 

24. 

n-Butyric acid 

<8 

35 

29.5 

26.4 


25. 

Pyridine 

spreads 

— 

— 


— 

26. 

ie 

Methylene iodide 

17.5 

75 

48 

43 

37.2 

27. 

Aniline 

25 

— - 

— 

— 

— 

28. 

* * 

Bisphenol A 

50 

— 

, 



62 . 2 


* 

Literature Value =50.4 

Melted at IGO^C and cooled to room temperature 
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TABLE-IV-C-6 
SUMMARY OF RESULTS 


Solid 


PTFE 


Surface Energy of 
Solids, erg/cm^ 


Literature Values 
erg/citi^ 


Dispersive 

Component 

Y^ = 19.6 

18.5^; 19.5^; 

19.6' 


Y^ =21.2 

23.9^; 22.6^; 

22.5' 



19.1^ 


Polar 

Component 

Y^ = 0.015 

2.0^ 


Y Total 

Y^ = 21.2 ± .35 

21.65 



HOPG 


Dispersive 

Component 


Polar 

Component 

Y Total 


33.5 

43.88 

41.7^ 


= 4.78 


Y^ = 48.6 ± 3.2 


33 ± 1.2 


8.1 ±3 
41.1 ± 3 


a,b 


Critical surface energy or surface tension 


dynes 


cm 


y, 


(6) (7) 


^'^Frora equation of state, Cos 0 = 2 (y „ /y )^ + 1^^®^ 

d 7 1 n \ 

Obtained from harmonic mean equation 

^Obtained from geometric mean equation 

^From liquid homologous, molecular weight dependence 

^Kaeble^^^ 

V fg) 

Drzal, L. T. , for HMU, i.e., high modulus fiber "un- 
surf ace- treated" 
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c). Discussion 

The driving force for the adhesion of a liquid on any 
substrate is the work of adhesion, i.e. , the reversible work 
required to separate a unit area of the liquid from the sub- 
strate. For the drop Of any liquid on HOPG (see Figure IV- 
C-11): 





Figure IV-C-11 

Contact Angle and Droplet Configuration 


where 




Y 


SL 


Yc 


Y< 


0 



= liquid-air surface tension 
= solid-liquid interfacial surface tension 
= solid-air surface tension 

= surface tension of HOPG in vacuum 
= the equilibrium spreading pressure 
= contact angle between liquid and HOPG 
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For air-HOPG equilibrium: 


0 


( 2 ) 


Using Young's equation, for HOPG-liquid-air equilibrium: 


Yj^Cos 6 = Ys - Ysl 


(3) 


Substituting (3) in (2) we obtain the Dupree-Young relation 
for work of adhesion: 

W = Yt (1 + Cos 9) + IT (4) 

a 'L e 

Assuming that as 6 ->■ 0, Yc 0 it 0 in orderly fashion, 

O 0 


W = Yt (1 + Cos 0 ) 
a. ij 


(5) 


The work of cohesion , i.e. , the work required to pro- 
duce two units of interfacial area from an original column of 
any liquid in air is given by 

= 2Yj, 

The difference between work of adhesion of a liquid on a sub- 
strate and the work of cohesion is the spreading coefficient: 

”a - ”o = - 2 Yi, = Sj^/s 

a) When W > W , S,. is positive, 9 = 0 and the liquid 

a c ±j/ b 

spreads spontaneously over the substrate to form a thin film. 

b) When W^ ' ^L/S coefficient) is negative, 

the liquid does not spread over the substrate but forms drop- 
lets or lenses with finite contact angle 9. 

c) When W = W , 

c a ' 

0 + 0 , Cos 9 = 1, =0 



14b 


d) When 0 = 180, 


«a = ®S/L = 


inefficiency factor = 



X 100 = 100. 


Table IV-C-7 lists values calculated using the expressions for 
work of adhesion, W , spreading coefficient, W - W and in- 

3. 3 C 

efficiency factor 


W - W 


2y. 


100 , 


as given above. 


Close observation of the values in Table IV-C-7 shows 
the work of adhesion of the liquid resins on PTFE is in the 
range of 46 to 52 erg/cm^ with a mean value of 48.1 erg/cm^ 
while for HOPG, W is in the range of 72 to 98 erg/cm^ with a 

3 

mean value of 83.7 erg/cm^ . The inefficiency factor may be 
considered as the fraction of the cohesive strength of the 
liquid that did not contribute to adhesion. These results 
strongly demonstrate that adhesional wetting on HOPG is not 
only vastly superior to that on PTFE, but that the work of 
adhesion can be enhanced by a high value of a high value 

of Yt f a low value for Yc-t » and a low value for spreading 

Jj bJ-i 

coefficient. Low spreading coefficient will occur when Yg 
is roughly about the same as Yt • 

Jj 

The critical surface energy y^^ of a solid is operational- 
ly defined by Zisman as 


Cos 9 = 1 + b(Y ” Yt ) 
' c L 


(6) 
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TABLE IV-C-7 

WORK OF ADHESION BETWEEN LIQUIDS, HOPG AND 
PTFE BASED ON CONTACT ANGLE MEASUREMENT 


Vlork of Adhesion 
Yt (1 + Cos 9) 

JLi 


Material [erg/cm^] 



HOPG 

PTFE 

Water 

109 

45 

Glycerol 

119 

58 

Epon 828 

83 

46 

Epon 871 

72 

50 

Epon 152 

78 

— 

Ciba MY720 

90 

48 

Ciba CY179 

79 

48 

Ciba 508 

£72 

52 

Tor Ion 4000 

98 

50 

Kelpoxy-G- 

272-100 

90 

47 


Spreading Inefficiency 


Coefficient 

S^ /o = w - W 
L/S a c 

[erg/cm^] 

Factoi 

W - W 
a c 

100 

* 

HOPG 

PTFE 

HOPG 

PTFE 

-36.5 

-100.5 

25 

69 

-26.5 

87.5 

21 

62 

-2.1 

39 

2.5 

46 

-0.5 

23 

0.7 

31 

-3.8 

— 

4.6 

- 

— 

-9.1 

-51 

9.2 

51.5 

-0.7 

-32 

0.9 

40 

— 

-26 


33.5 

-21.6 

-69 

18 

57.5 

-7.4 

50.4 

7.5 

51.8 


* 

Y 


= Mean Yt predicted from PTFE and HOPG 


L 
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v/here b is a constant, and Yt is the surface energy of the 
liquid. In all instances in which Yjv < f the equation pre- 
dicts the contact angle will be zero,' and the liquid will 
spread. Also, that for all y^ ^ Y^ > the cosine of 6 decreases 
in a linear manner. 

Substituting (6) in (5) , the work of adhesion becomes 

w; = (2 + bY^)YL - bY^ 


which is a parabola. Equation (7). predicts that w' will dis- 

1 • T « 4. dW' ,, 1 \ 1 

play a maximum value at = 0, Y^ = g + 2 '^ c 

L 

and that the maximum work that can be obtained from this sur- 


face will be given by 

w;(max) =1+1 + jbY^ (8) 

A plot of maximum work obtainable from an HOPG surface with 
Y^ , critical surface tension, for various values of b, shows 
that for W(Y|^ = 33.5) = 80.3 erg/cm^ and W(y^ = 48.6) - 103,2 
erg/cm^ for b = 0.0252, the work factor improved by 28 percent 
for a 45 percent increase in y^ (Figure IV-C-12) . This range 
compares favorably with the values in column 2 of Table IV-C-7 
for the various epoxies. Equation (8) suggests that the only 
ways more work can be obtained from the HOPG basal plane sur- 
face are 

a) decreasing the value of Zisman slope b and 

b) increasing the critical surface energy of 
HOPG by oxidative activation of HOPG surface 
so that bonding reactions can occur with 
resin, etching or other techniques that may 
raise y^ 01 ^ the interaction between liquid 
and solid . 



W(yc) = 


b= 1/400. 


b= 1/320 


b= 1/240 


b = l/l60 


160 


b=l/4a 




zj Epoxy Resins 


O 20 40 60 80 100 120 

erg /cm ^ ► 

Figure IV-C-12 

Plot of tlaxiiauin Work W (max) Obtainable on HOPG at 
Various Critical Surface Energy for Different b Values 
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In fact, it could be argued that b and related and 

that any factor that will affect b also affects The dif- 

ficulty with etching HOPG surface or HM fibers to improve y 

c 

is that some of the highly oriented outermost layers are re- 
moved and also additional flaws may be introduced on the sur- 
face of the fibers. Hence, shear strength will increase at 
the expense of tensile strength of the composite. 

The scatter y =20 erg/cm^ in the Zisman plot could be 
traced to the following assumptions and reasons: 

a) Liquids identical in chemical properties to 
a solid would have a zero contact angle on 
the solid, hence the same surface tension. 



b) All liquids with identical y^^ will display 
the same contact angle on any solid (con- 
verse of assumption "a"). 

c) Any liquid whose surface tension is below 
the critical surface tension of a solid 
will completely wet the solid. 

d) There is a varying degree of adsorption of 
the various liquids on the HOPG surface, 
which leads to a hysteresis. This could be 
shown vividly in the case of diethylene- 
triamine on HOPG. (Figure IV-C-13) 

e) Each surface of HOPG is slightly different, 
in terms of roughness of basal planes and 
'fractional coverage' of non-basal planes 
with edges or high surface energy regions 
that could not be observed at 10 OX. 
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Figure XV-C-13 

Drop of Diethylenetriamine on HOPG Surface 
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f) The polar surface energy components of the 
various liquids and HOPG were ignored. 

The Zisman relation works well for homologous series of 

liquids, such as alcohol as shown in Figure IV-C-10. Hov;- 

ever, assumptions a) and b) above become particularly poor 

for non-homologous liquids. This is because it is the 

matching between both the polar and dispersive components of 

the liquid and HOPG that determines the contact angle. The 

empirical method of estimating may result in a value that 

d 

may agree well with rather than the true solid surface 
tension y^ , particularly where substantial contribution from 
y^ exists. 

O, 

On non-polar and low energy surfaces the third assump- 
tion, c) , suggests that complete wetting can only occur when 
the surface tension of the wetting liquid has been reduced to 
a critical value, y^ , characteristic of the substrate, e.g. , 
33.6 dynes/cm^ for HOPG and about 31 dynes/cm^ for polyethyl- 
ene. This equates y with y„ - yet • For many liquids, a 
plot of Cos 9 versus y^. is linear. Extrapolation of this 
plot to Cos 9=1, that is 9 = 0, yields y^ or the critical 
surface tension for complete wetting. This quantity is a 
constant for a particular substrate irrespective of the nature 
of the wetting agent. This will, of course, require that 

(8) 

Yo “ Yot reduced to the same value in all cases. Rosen' ' 
o oJLi 

has suggested that the two necessary conditions to achieve a 
constant value of yg - Ygj^: (a) the free energy per unit 
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area,Yc^ of surface in equilibrium with the liquid- 

saturated air above it be equated with Yg / ths free energy 

of the solid, free of adsorbed material (This implies that 

TT^ , the equilibrium pressure, will have to vanish.) ; and (b) 

the adsorption of liquid on the substrate interface occurs 

with some orientation and degree of packing, producing the 

same value Yot (thus a constant value of y„ “ Yc?t ) the 

surface tension of the spreading liquid approaches y. 

c 

Schwarz in 1964 and Bernett^^^^ in 1959 have shown 

that the presence of highly fluorinated carbocyclic acids and 

their salts will lower the value of y for polyethylene from 

c 

its usual value of almost 31 dynes/cm^ to about 20 dynes/cm^ . 
This adsorption of the fluorinated surfactants onto the poly- 
ethylene surface gives the result that solutions of these 
surfactants having surface tensions less than the normal y^ 
for polyethylene do not spread on it. 

In our work, the contact angle of Epon 828 on HOPG is 
about fifteen degrees, while that of DTA is thirty degrees. 
One would have expected some intermediate contact angle in 
the mixture of these two materials. In fact, the contact 
angle remains consistently at thirty degrees for all mixtures 
because the amine preferential migrates to the interface and 
spreads out. Although the contact angle of distilled water 
on HOPG is sixty degrees, in the case of a mixture of Epon 
828 and distilled water, the contact angle of the mixture re- 
mains at fifteen degrees because the water phase separated as 
small bubbles in the droplet. 
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The requirement that the surface tension of the wetting 
liquid be reduced by the surfactant to some critical value 
characteristic of the substrate is thus necessary but not 
sufficient for complete spreading-wetting. In general, a 
liquid will not produce complete wetting of HOPG if its sur- 
face tension is above the critical surface tension, y . For 
liquids with surface tension below 33.6 dynes/cm^ , complete 
wetting of HOPG may or may not be produced as in the case of 
n-butyric acid, acetic anhydride and acetonitrile. 

In the fourth assumption, (d) , the effect of surface 
roughness could be relaxed provided the plane remains basal; 
for instance, no significant contact angle change was observed 
on growth cones on HOPG (basal) . The effect of non-basal 
planes was observed to be quite drastic. Contact angle has 
been observed to be reduced by as much as nineteen degrees 
by the presence of edges. The liquid tended to be dragged 
along this edge. 

d) Conclusion 

Our preliminary work in the study of surface energetics 
of the basal planes of HOPG suggests that the polar contribu- 
tion to the surface is as much as ten percent. The critical 
surface energy of HOPG, y^ , and Zisman slope, b, have been 
characterized, and we have shown that the calculated work of 
adhesion for various epoxies on HOPG is quite close to the 
maximum work obtainable from the surface. 
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3. Plans for the Upcoming Period 

While the separation of the surface tension of HOPG into 
its polar and dispersive component indicates that at 20-25°C 
the polar contribution accounts for roughly ten percent of 
the total surface energy, in HMU fiber the polar contribution 
accounts for as much as twenty percent, according to the work 

(9) 

of Drzal et al. This separation of the surface energy 
component for the various epoxy and resin matrices would re- 
quire about five different solid surfaces with characterized 
polar and dispersive components, respectively. This is one 
of our immediate objectives. 

The relationship between contact angle of resins and 
their viscosity is not well know at this time. The separa- 
tion of surface tension of epoxy into their polar and dis- 
persive components, the relationship between contact angle 
and viscosity of resins and the various surface modifica- 
tions to lower b and to increase the critical surface ten- 
sion is to be performed. 
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IV-D FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator: E. Krempl 

The deformation and fatigue failure behavior of compos- 
ites is being studied for (.essentially one-dimensional) lami- 
nate samples under uniaxial loading and, more intensively, 
for graphite-epoxy tubes under biaxial (tension, torsion) 
loading. The aim of this research is to provide basic under- 
standing and design information on the biaxial response of 
advanced composites. 

1. Status 

The following had been accomplished prior to the start 
of the current reporting period: a series of fatigue tests 

to failure were performed, using uniaxially loaded laminate 
specimens, to evaluate frequency dependence and the influence 
of R - ratio. Specially designed test fixtures were manu- 
factured capable of loading tubular specimens biaxially and 
triaxially (axial, torsional and internal pressure loading). 
Eight tubular graphite-epoxy specimens with four layers of 
(±45°) were designed and manufactured. The elastic proper- 
ties of these specimens were determined, using strain measure- 
ment on the gage length. 

Studies of damage accumulation models have also been 
performed with particular reference to the residual strength 
degradation model. Models in the literature have been compared 
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for such aspects as constant amplitude loading, rate depend- 
ence and non-deterministic properties, etc. 

2. Progress During Report Period 

Additional tubular test specimens have been manufactured, 
and their elastic properties under uniaxial , torsional and 
combined loading have been determined. Elastic properties 
were determined using an Instron biaxial extensometer attached 
to the gage length of the specimen under cyclic axial, tor- 
sion and combined axial-torsion loading, covering the first 
and third quadrant of the tension-torsion space. 

Figures IV-D-la and-lb show that the slopes of axial 
load versus axial strain and torque vs. engineering shear strain 
are unaffected by proportional loading. In the elastic range 
the material behaves isotropically. Figures IV-D-la and -lb 
also demonstrate that the moduli are equal in tension and 
compression. 

The small hysteresis loops in Figures IV-D-la and -lb 
are probably due to friction in the extensometer. They do 
not appear to be a property of the tubes. 

With all tubes we measured effective Poisson's ratio 
using an axial and a diametral MTS extensometer. The effec- 
tive Poisson's ratio (negative diametral strain/axial strain), 
the effective Young's modulus and the effective shear modulus 
are listed in Table IV-D-1. 

The post-elastic behavior of these tubular specimens has 
also been determined, including ultimate strength under 
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Figure IV-D-lb 
(companion to Figure -la) 

Shear Stress Vs. Engineering Shear Strain (Angle Change) 
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TABLE IV-D-1 



STATIC PROPERTIES OF 

GRAPHITE -EPOXY TUBES 

(±45°) 

s 

Tube 

No. 

Modulus 

E (GPa) G (GPa) 

Ultimate 
0u (MPa) 

strength 
TU (MPa) 

Poisson 

Ratio, 

7 

12.92 




' — _ 

9 

11.84 

27.20 

* 

-148 



0.784 

10 

12.00 

22.74 

* 

-166 

. — — — 

0.736 

11 

12.93 

23.56 

** 

148 

— 

0.724 

12 

11.40 

23.14 


18 0 

0.752 

13 

12.31 

23.41 



0.738 

14 

13.83 

27.44 



0.748 

15 

11.34 

___ 



0.875 

16 

11.76 




0.722 

17 

11.34 




0.884 

18 

11.34 




0.696 

20 

11.13 




0.863 

21 

12.86 




0.447 

22 

12.00 




0.769 

23 

11.88 




0.745 

Avg. 

11.99 

24.58 



0.749 

* 

Compression 





** 







Tension 
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tension, compression and torsional loading. 

Figures IV-D-2a and -2b show graphs of axial load versus 
crosshead displacement. Torque versus angle of twist of the 
actuator is plotted in Figure IV-D-3. For axial loading, the 
diagrams depart from linearity at 26 percent and at 35 percent 
of the ultimate in tension and compression, respectively. The 
nonlinear portions are characterized by time-dependent be- 
havior as evidenced by the increase of the displacement (creep) 
during periods of constant load (see Figure IV-D-2a) . 

Fatigue testing has also been conducted, using completely 
reversed, load-controlled uniaxial loading. Axial and diamet- 
rical strain measurements were taken during these tests. Fig- 
ure IV-D-4 shows behavior during such loading at a load ampli- 
tude of about 60 percent of the ultimate load. This demon- 
strates further the time dependence shown in Figure IV-D-2a. 
Cycles 9 and 10 (ABCD) took 240 seconds. When B was reached 
during the eleventh cycle the rate was reduced by a factor of 
1000. A "creeping motion" , BE , is evident before the load in- 
creases again. The strain range of the "slow" loop (BEFGH) , 
which takes 66.67 hours for completion, is considerably larger 
than that of the "fast" loop (ABCD) . The strain range in- 
creased further during the next loop which started at H and 
ended at I. 

The torque versus angle of twist of the actuator curve. 
Figure IV-D-4, shows much less deviation from linearity than 
the axial curve due to the "fiber dominance" of this defor- 
mation. Deviations from linearity are beyond 50 percent of 
the ultimate load. 
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As indicated on Figures IV-D-2a, -2b and -3, the ultimate 
load was determined in each case. The appropriate ultimate 
strengths are listed in Table IV-D-1. 

It is noteworthy that the tubes permit the determination 
of the compressive ultimate strength without any buckling 
problems. The ultimate compressive strength would appear to 
be very comparable to the ultimate tensile strength. 

3. Plans for the Upcoming Period 

The completely reversed, load-controlled axial tests 
which are now under way to ascertain the influence of loading 
rate on the hysteresis loop and on the fatigue life will be 
continued and completed. They are to be followed by torsion 
and combined axial-torsion tests. 

A paper on biaxial deformation and fatigue behavior is 
scheduled to be written during the next reporting period. 
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IV-E ACOUSTIC EMISSION OF COMPOSITE MATERIALS UNDERGOING 
FATIGUE 

Senior Investigators: H. Otsuka 

H. A. Scarton 


The Acoustic Emission (AE) technique has been shown to 

be useful in monitoring composite material fracture under 

(1 2 ) * 

conditions of monotonically increasing loads . Differ- 

ent signal amplitudes were shown to be associated with dif- 
ferent AE sources; namely, a) low amplitude signals, matrix 
crazing; b) intermediate amplitude signals, delamination and 
fiber pull-out and c) large amplitude signals, fiber break- 
age. The current research deals with the extension of this 
work to include fatigue failure. 


1. Status 

Acoustic emission test and data analysis have been per- 
formed for composite material tensile specimens under con- 
stant load rates. Both count versus time data and length- 
location of acoustic emission sources have been successfully 
obtained for (±45°) and (0° -90°) graphite cloth epoxy compos- 
ites. 

2. Progress During Report Period 

The shape of the specimens used in fatigue tests is 
shown in Figure IV-E-1. Reinforcing aluminum tabs were bonded 

* 

Niimbers in parentheses in this section refer to the refer- 
ences listed on page 175. 
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Figure IV-E-1 


Test Specimen (dimensions in inches) 
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to the composite specimens as shown. Load was transferred 
from the testing machine to the specimens through mechanical, 

85 dB attenuator fixtures by steel pins. The materials used 
in the specimens were carbon cloth-ply with a 0°-90° fiber 
orientation and a twelve-layer prepreg, carbon-epoxy compos- 
ite with the following stacking sequence: (0°, ±45°«, 90°) . 

Testing was performed on an Instron 1333 (capacity 50 kN, 

25 kN range) using the AE setup shown in Figure IV-E-2 . 

Amplitude distributions of the acoustic emission signals 
were recorded in memory A and location distributions in mem- 
ory B. 

Tension-tension fatigue tests were performed using both the 
0°-90° and twelve-layered graphite-epoxy composites according to 
the load schedule shown in Figure IV-E-3 , where the specific 
loads undergone and fatigue strength are indicated in Table 
IV-E-1. Figures IV-E-4 and -5 show the AE count versus fa- 
tigue cycles for the 0°-90° and twelve-layer prepreg speci- 

( 2 ) 

mens, respectively. The previously reported static tests 
showed sharply increasing AE signals (critically active) just 
prior to failure. The fatigue results presented here show 
(after an initial period) a linear relationship between fa- 
tigue cycles and AE count. This suggests that differences in 
failure mechanism may exist between these two forms of fail- 
ure. Critically active behavior was observed on the 0°-90° 
specimens in fatigue just prior to failure, while no such 
activity was observed on the prepreg specimen. No defect was 
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Figure IV-E-2 

Acoustic Emission Instrumentation 
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TABLE IV-E-1 
FATIGUE TEST RESULTS 


Load 

Number of Mean Value ± Variation 

Specimen Cycles (psi) (kg/mm^ ) (% of u.T.S.) 


#6, 0-90“ 










0->300 

43390 

± 

5784 

30.27 

± 

4.034 

75.00 

+ 

10 

300^600 

44130 

± 

5784 

30.78 

± 

4.034 

76.28 

± 

10 

600^900 

44860 

± 

5784 

31.29 

+ 

4.034 

77.55 

± 

10 

900^1200 

45610 

± 

5784 

31.81 

± 

4.034 

78.83 

± 

10 

1200->1500 

46340 

± 

5784 

32.32 

± 

4.034 

80.10 

± 

10 

1500^1800 

47080 

+ 

5784 

23.84 

± 

4.034 

81.38 

+ 

10 

1800->2100 

47810 

± 

5784 

33.35 

± 

4.034 

82.65 

± 

10 

2100^2400 

48560 

± 

5784 

33.89 

± 

4.034 

83.93 

+ 

10 

2400^2700 

49290 

± 

5784 

34.38 

± 

4.034 

85.20 

± 

10 

2700-^3000 

50030 

± 

5784 

34.89 

± 

4.034 

86.48 

± 

10 

3000^3300 

50770 

± 

5784 

35.41 

± 

4.034 

87.75 

± 

10 

3300^3600 

51500 

± 

5784 

35.92 

± 

4.034 

89.03 

± 

10 

3600^3862 

52250 

± 

5784 

36.44 

± 

4 . 034 

90.30 

± 

10 


#9, 0-90“ 



0 

/ The 

same as 

for Specimen 

. #6. 





3600 

3900^3997 

52980 

+ 

5784 

,36.95 

± 

4.034 

91.58 

.+ 

10 

P-1 

12 Layer 

0^29 

49194 

± 

9838 

•34.31 

± 

6.86 

83.33 

± 

16.67 

P-2 

12 Layer 

0^300 

41250 

± 

8250 

28.77 

± 

5.940 

69.88 

± 

13 .98 

300->600 

44000 

+ 

8250 

30,69 

± 

5.940 

74.54 

± 

13.98 


600->-669 

49500 

± 

8250 

34.53 

± 

5.940 

83.85 

± 

13 . 98 

P-3 

12 Layer 

0-^300 

41250 

± 

8250 

28.77 

± 

5.940 

69.88 

+ 

13.98 

300->322 

46750 

+ 

8250 

32.61 

± 

5.940 

79.19 

± 

13.98 

P-4 

12 Layer 

O-J-150 

41250 

± 

8250 

28.77 

± 

5.940 

69.88 

± 

13 . 98 


150^300 

42350 

± 

8250 

29.54 

± 

5.940 

71.74 

+ 

13.98 


300^450 

43450 

± 

8250 

30.31 

± 

5.940 

73.60 

± 

13.98 


450^600 

44550 

± 

8250 

31.07 

± 

5.940 

75.47 

± 

13.98 


600^750 

45650 

± 

8250 

31.84 

+ 

5.940 

77.33 

± 

13.98 


750-^900 

46750 

± 

8250 

32.61 

+ 

5.940 

79.19 

± 

13.98 


900-»-1036 

47850 

± 

8250 

33.38 

± 

5.940 

81.06 

+ 

13.98 

P-5 

12 Layer 

0^150 

41250 

± 

8250 

28.77 

± 

5.940 

69.88 

+ 

13.98 


150^300 

42350 

± 

8250 

29 . 54 

± 

5.940 

71.74 

± 

13.98 


300-5-450 

43450 

± 

8250 

30.31 

± 

5.940 

73.60 

± 

13.98 


450-5-600 

44550 

± 

8250 

31.07 

± 

5.940 

75.47 

± 

13.98 


600^750 

45650 

± 

8250 

31.84 


5.940 

77.33 

+ 

13.98 


750-5-900 

46750 

± 

8250 

32.61 

± 

5.940 

79.19 

± 

13.98 


900-5-990 

47850 

± 

8250 

33.38 

± 

5.940 

81.06 

+ 

13.98 



AE-TOTAL COUNT 





AE-TOTAL CCXINT 
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Figure IV-E-5 

Acoustic Emission Counts versus Load Cycles 
for the Twelve-Layer Carbon Prepr eg Specimen 
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apparent on the 0°-90® specimen, while the twelve-layer spe- 
cimen failed as a result of the growth of an easily visible 
defect at the center 0°-90° layers (Figure IV-E-6) . The 
ultimate tensile strengths of the 0® -90 “and twelve-layer 
specimens were 57,851 psi (0.399 GPa). and 59,032 psi (0.407 
GPa) , respectively. Amplitude distributions (Figure IV-E-7) 
and location distributions (Figure IV-E-8) were recorded 
every 150 cycles during the fatigue test. It is possible 
that the transitions from one given amplitude distribution to 
another correlate to the growth of the 0® -90 "center layer 
defect. Note that the fixed nature of the multi-modal distri- 
bution is maintained with increasing load. 

3 . Plans for the Upcoming Period 

Further work will be devoted to analyzing the acoustic 
emission process as it relates to the fatigue process as a 
method of evaluating composite material performance. 

4. References 

1) Otsuka, H. and H. A. Scarton, "Variation in Acoustic Emis- 
sion between Graphite- and Glass-Epoxy Composites", Journal 
of Composite Materials, Submitted for publication July 1980 

2 \ Otsuka, H. and H. A. Scarton, "Acoustic Emission Testing of 
Composite Tensile Specimens", NASA/AFOSR 38th Semi-Annual 
Report, "Composite Structural Materials", June 1980, pp 133- 
139. 
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Figure IV-E-7 

Transition of Amplitude Dis- 
tribution for Twelve-Layer 
Specimen (#5) with Number of 
Cycles 
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Figure IV-E-8 

Transition of Location Distribution for Twelve-Layer Specimen (#5) with Number of Cycles 
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IV-F RESIN MATRIX DEPENDENT CHARACTERISTICS OF COMPOSITES 

Senior Investigator: S. S. Sternstein 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
physical properties of the matrix resin. To date, the re- 
search has involved the precise viscoelastic characterization 
of epoxy neat resins, interlaminar failure of composites and 
the inhomogeneous swelling of and the effects of moisture on 
composites. 

1 . Status 

Dynamic mechanical spectroscopy, creep and stress relax- 
ation techniques have been used to determine the effect of 
moisture and temperature on the properties of neat resins. 
Non-linear plate theories have been reviewed to guide planning 
of tests intended to determine interlaminar failure tolerance 
of the epoxy layers between adjacent plies. Constitutive re- 
lations have been formulated describing the thermodynamics 
of water interaction with epoxy resins. The model includes 
a water dependent modulus of elasticity. Solutions for 
stress and strain around a carbon fiber in an infinite epoxy 
matrix have been obtained. 

2. Progress During Report Period 

Our activities during the current period have emphasized 
three areas, namely, a) the delcimination failure of carbon- 
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epoxy laminates as a result of out-of-plane loading; b) the 
stress distribution in various stacking sequence laminates 
subjected to biaxial, out-of -plane bending; c) the viscoelas- 
tic characterization of situ resin as a function of both 
humidity and temperature. A fourth project on resin cure 
rheology has been initiated with non-NASA support and is also 
relevant to this program. 

a) Delamination Failure 

This project has emphasized the development of a simple 
but effective test for the study of delamination failure in 
carbon-epoxy laminates. The test sample consists of a circu- 
lar disc of roughly 1.3 inches diameter and approximately 10 
to 14 plies thickness. The sample is mounted on a special 
anvil which is attached to an Instron compression load cell. 
The anvil supports the disc sample around the perimeter at a 
one inch diameter. The out-of-plane load is applied to the 
center of the sample using a large ball bearing which is 
loaded by the Instron movable crOsshead. The test geometry 
may be described as a centre symmetric bend test. It is im- 
portant to precisely center the ball bearing with respect to 
the sample anvil so as to maintain a true biaxial deformation 

The load cell circuitry of the Instron ha& been modified 
to permit external recording of the load signal. This neces- 
sitated building an amplifier, filter and buffer circuit. 
Crosshead motion is monitored using a highly sensitive LVDT 
and external amplifier. Thus, direct x-y recording of force 



181 


and displacement can be obtained with high accuracy. Finally, 
the compliance of the machine and load cell have been care- 
fully measured and used to correct measured sample compliances 

The centre symmetric bend test is useful for delamination 
failure measurements and has been shown to give very reproduc- 
ible data. It is important to note that this test always 
initiates failure at the center of the sample and not at the 
edges. Therefore, the presence of nicks and flaws due to 
sample preparation is of no significance. While it has these 
important advantages, however, the centro symmetric bend test 
should not be viewed as providing a polar-symmetric situation. 
The fact that each filament layer has non-zero thickness 
means that the three important moments of lamina thickness 
about the plate neutral axis (zero, first and second order of 
the depth dimension) cannot all be made identical, as one 
changes stacking angle from one lamina to another . 

Samples were prepared in a press from prepregs. The 
stacking sequences used thus far include the following 12-ply 
laminates: (A) [0/90]2g/ (B) [0 , 45 , 90 , 135 , 0 , 45] g , and (C) 

I0,30,50,90,120,150]g. In choosing these, priority was given 
to perserving the angle between adjacent plies through the 
depth of a particular sample. Sample A will be referred to 
as a 90° sample in that each ply rotates by 90°; similarly, 
samples B and C are 45° and 30° samples, respectively. 

The stacking sequences all contain a screw axis symmetry 
of 90°, 45° or 30° per ply. Note that sample B, which has 
more layers with filaments in the +45° direction than the -45° 
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direction could only have been "balanced", in that sense, by 
either relaxing the "equal angle between adjacent plies cri- 
terion" or changing the total number of plies from that used 
for samples A and C. Retaining the relationships among sam- 
ples was considered more important than "balance" in that 
conventional sense. 

A typical load versus disc (center) deflection curve is 
shown in Figure IV-F-1 for a 90® sample. The characteristic 
vertical drop around 80 kg load is the first delamination 
failure. Fresh samples either unloaded or loaded to less 
than roughly 70 kg have been examined in the Acoustic Imaging 
Lab of Professor Das at R.P.I. and show no signs of delamina- 
tion; conversely, the sample shown in Figure IV-F-1 clearly 
shows a large central delamination when examined ultrasonic- 
ally. 

Figure IV-F-2 shows another 90° sample which was first 
subjected to 50 cycles of deflection between 0 kg and 50 kg. 
The sample was then tested to 100 kg. Clearly, the precycle 
test had little effect, if any, on the 80 kg load drop. 

A 90° sample was preconditioned in boiling water for 
1000 hours and then tested, as shown in Figure IV-F-3 . The 
load at which failure occurs is markedly reduced to about 
35 kg, but the sharp load drop as seen in Figure IV-F-1 is 
now an inflection with a significant drop in stiffness be- 
yond that point . 

The first cycle failure of a dry 45° sample is shown in 
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Figure IV-F-1 

Lioad vs. Deflection at Disc 
Denter (90° sample, no pre- 
conditioning) 
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Figure iy-F-4 . When compared to Figure IV-F-1 , it is seen 
that the delamination failure point is reduced to roughly 
67 kg (ys% 80 kg for the 90° sample) and the load drop is 
much smaller. 

A brief summary of a statistically significant niimber of 
samples is given in Table IV-F-1. 


TABLE IV-F-1 

FAILURE CHARACTERISTICS OF VARIOUS LAMINATES 


Property 

PgCkg) 

^g(kg) 

Xg (mm) 


90° 

Sample 

79.7 

15.0 

0.683 


45° 

Sample 

68.6 

2.70 

0.619 


30° 

Sample 

70.4 


Inflection 

0.669 


where: P_ = mean load at first failure 

^ = mean load drop at first failure 

13 

X„ = mean deflection at first failure 


Finally, we note that the results are- reasonably inde- 
pendent of sample size. Figure IV-F-5 shows a 90° sample 
which was tested using an anvil /2 times the diameter used in 
Figure IV-F-1. The first delamination failure occurs at 87 kg 
and has a load drop of 14 kg which is within the statistical 
expectation for 90° samples. 


b) Stress Analysis 

A detailed analysis of the elastic stress distribution 
in the laminates used in section a) above is being performed 
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Figure IV- 

Load vs. Deflectic 
Center (45° sample 
ditioning) 






W, (mm) 


Figure IV-F-5 

Load vs. Deflection at Disc 
Center (90® sample, no precon- 
ditioning; outer dia. times 
standard) 
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using modifications of laminated plate theory. In particular, 
we have been concerned with the state of stress at the inter- 
faces between plies. For this reason it has been essential 
to satisfy both the condition of stress equilibrium at each 
interface as well as the condition of displacement continu- 
ity. Simple laminated plate theory as used by many authors 
does not give equilibrium of the stresses at each interface. 

To provide such equilibrium, at least a quadratic distribu- 
tion of displacement through each ply thickness is called 
for; a different quadratic is required for each ply. 

Some results are shown here, but a complete analysis and 
discription of our findings will be published shortly. The 
out-of-plane stress component a (rz) is shown versus thickness 
(depth z) for a 90° sample in Figure IV-F-6 at angle 0°. The 
same stress component, cr (rz) , is shown at angle '30° in Figure 
IV-F-7 . Similarly, the stress component, a(c))z), is shown at 
an angle of 30° in Figure IV-F-8 . 

The maximum shear stress lout-of -plane resultant of 
and 0 {rz)J is shown versus angle for the various interfaces in 
Figure IV-F-9 and a polar representation is given in Figure 
IV-F-10. 

The maximxam out-of -plane shear stress resultant for a 
sample having a stacking sequence of (0)^2 shown in Figure 
IV-F-11 and should be compared to Figure IV-F-9 for the 90° 
sample . 
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Figure IV-F-10 

Polar Representation of Maximum Shear Stresses for Various Interfaces Between Lamina 

(90° sample) 
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Figure IV-F-11 

Maximum Shear Stress vs. i^gle for Various Interfaces Between Lamina 
(unidirectional, sample) 
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c) Viscoelastic Study 

The sample and loading geometry described in section 
a) has also been used in our dynamic mechanical spectrometer 
to study the situ glass transition behavior of the epoxy 
resin. The sample is under load servo control and is pro- 
grammed with a static load of 3 kg and a dynamic load of 1 kg 
amplitude. Frequency and temperature are computerized and 
the dynamic analysis (dynamic modulus and loss factor) is 
performed using a Fourier type correlation in an instrviment 
referred to as a "Dynalyzer". 

Typical results ..for a 90° sample are given in Figures 
IV-F-12 and 13. Figure IV-F-12 shows the loss factor (tan 
delta) at 6.5 Hz versus temperature. The effect of water is 
quite apparent after the sample is reheated. Notice the re- 
duction in mechanical loss at 100°C and increase by 5° of 
the glass transition (at 160°C) upon the second heating. 

This is not due to post-curing. The dynamic modulus Ie*| is 
shown in Figure IV-F-13 . These "small-strain" dynamic re- 
sults are extremely useful in the characterization of samples, 
particularly in regard to sample batch-to-batch variations in 
either processing or starting materials. 

3. Plans for the Upcoming Period 

We intend to continue the centro symmetric bending tests 
using additional ply arrangements to further ascertain the 
effects of temperature, humidity and load rate. Further, an 
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Figure IV-F-13 

Dynamic Modulus |e*| vs. Temperature (®C) 
(90° sample at 6.5 Hz) 
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attempt will be made to extend this testing technique to the 
evaluation of neat resins. The dynamic mechanical spectros- 
copy results are being obtained at frequencies from 0.1 Hz to 
40 Hz and master curves will be constructed to ascertain the 
applicability of time-temperature superposition to the in 
situ glass transition behavior of the epoxy and also the re- 
laxation time spectrum of situ epoxy relative to neat 


resin. 
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